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SECTION  I 


INTRODUCTION 

Because  of  its  unique  rotation- vibration  structure  and  its  domi- 
nant role  in  atmospheric  absorption  of  C02  laser  radiation,  the  H2160 
molecule  has  recently  been  a topic  of  great  interest  (refs.  1-10). 

The  internal  nuclear  motion  of  this  molecule  is  substantially  complicated 
by  the  following  four  features:  (1)  the  large  molecular  asymmetry; 

(2)  the  small  moments  of  inertia  (and  correspondingly  large  rotational 
energies);  (3)  the  extreme  anharmonic  and  low  frequency  structure  of 
the  bending  mode  vibrations*;  and  (4)  the  large  amplitude  fluctuations 
of  the  bond  angle.  These  features  cause  water  to  display  strong 
coupling  between  the  molecular  rotations  and  the  bending  mode  vibration 
and,  in  addition,  to  exhibit  large  centrifugal  distortion  effects 
(ref.  1).  These  effects  manifest  themselves  not  only  in  the  rotation- 
vibration  energy  level  scheme  of  water  but  also  in  the  effective  geo- 
metric structure.  For  example,  examination  of  recent  high-resolution 
microwave  data  (refs.  1,2)  reveals  that  the  apparent  bond  angle  deviates 
over  a range  of  3.4°  depending  on  which  pair  of  (experimentally  de- 
termined) inertia  constants  are  used  in  the  analysis.  (The  apparent 
bond  length  deviates  over  a range  of  0.014  K,  which  is  appreciably  less 
than  that  of  the  bond  angle.) 

Present  analysis  of  high-resolution  microwave  and  far-infrared 
spectral  data  involves  a least  squares  fitting  to  a phenomenological 
rotational  expansion  Hamiltonian,  first  developed  by  Watson  (ref.  11). 

For  virtually  all  molecular  systems  of  interest,  such  spectroscopic 
modeling  techniques  yield  a reliable  and  useful  description  of  the 
molecule.  However,  these  techniques,  which  rely  on  the  rapid  convergence 
of  a power  series,  require  that  the  rotation-vibration  interaction  and 
centrifugal  distortion  be  in  some  sense  small.  For  the  particular  case 

* The  experimental  energy  spacings  are  1594.5,  1556.94,  and  1515.25  cm"1 
(from  refs.  3 and  17). 
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of  water,  this  is  not  the  case.  Thus,  one  is  lead  to  question  the 
validity  of  spectroscopic  modeling  techniques  as  applied  to  water.  In 
this  report,  this  question  is  examined  in  detail. 

Before  describing  the  nature  of  this  work,  a brief  statement  is 
presented  on  the  results.  As  is  well  known,  the  standard  Watson  rota- 
tional Hamiltonian  consists  of  a series  expansion  in  the  components  of 
the  angular  momentum  operator  to  various  powers  together  with  a number 
of  rotation  and  rotation-distortion  coefficients  that  are  to  be  determined 
from  available  experimental  data.  For  light  asymmetric  molecules  such 
as  water,  the  rotational  wave  functions  depend  in  a sensitive  way  on 
these  coefficients.  Hence,  if  one  requires  accurate  wave  functions  (and 
therefore,  a correct  physical  description  of  the  molecule),  it  is  im- 
portant that  these  coefficients  be  unique  and  that  the  power  series 
expansion  be  highly  convergent.  These  questions  have  been  examined  for 
two  systems:  (1)  a rotating  Morse  oscillator  and  (2)  the  H2O  molecule 

itself.  We  find  that  for  the  high-lying  states  (i.e.,  J s 10)  the 
coefficients  are  not  unique  (i.e.,  they  vary  with  the  levels  used  for 
fitting),  the  series  expansion  does  not  converge,  and  the  higher  order 
coefficients  are  devoid  of  any  physical  significance.  Furthermore,  the 
resultant  rotational  Hamiltonian  has  little  value  with  regard  to  its 
ability  to  either  predict  new  energy  levels  (i.e.,  states  that  were  not 
used  for  fitting  purposes)  or  to  determine  rotational  wave  functions 
that  accurately  describe  the  physical  state  of  the  molecule.  As  a 
consequence,  this  approach  appears  to  be  essentially  useless  for  describing 
the  high-lying  states  of  H2O.  These  difficulties  arise  from  and  reflect 
the  large  rotation- vibration  interaction  in  water  and  are  discussed  in 
detail  in  this  report. 

In  two  companion  papers,  which  were  guided  and  strongly  influenced 
by  the  findings  of  this  work,  the  rotation-vibration  structure  of  the 
water  mulecule  is  calculated  in  a way  that  directly  includes  the  bending 
mode. 
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SECTION  II 


GENERAL  DISCUSSION 

It  is  well  known  that  a theoretical  analysis  of  a molecular  sys- 
tem undergoing  both  rotational  and  vibrational  motion,  even  within  the 
Bom-Oppenheimer  approximation  (ref.  12),  is  usually  difficult.  This 
is  especially  true  for  moleculer  systems  that  exhibit  strong  nonlinear 
vibrational  features  or  significant  centrifugal  distortion  effects, 
i.e.,  light  molecules  such  as  water. 

A number  of  approaches  may  be  applied  to  the  study  of  the  water 
molecule.  The  first  would  be  to  determine  the  electronic  wave  functions 
and  to  use  these  to  calculate  the  vibrational  potential.  Using  this 
potential,  the  nuclear  wave  equation  would  then  be  solved  to  obtain  the 
rotation-vibration  energies  and  wave functions . Unfortunately,  the 
large  number  of  first-principle  calculations  (ref.  13)  of  the  molecular 
electronic  structure  of  water  have  not  produced  wavefunctions  of  suffi- 
cient accuracy  to  render  this  approach  useful. 

The  second  approach  is  to  choose  a nuclear  potential  that  will 
fit  the  experimentally  observed  vibrational  (e/  = 0 -**  J = 0)  energy 
differences.  This  potential  is  then  combined  with  the  Wilson-Howard 
(ref.  14)  expression  for  the  molecular  kinetic  energy,  and  the  resulting 
Hamiltonian  is  used  to  obtain  the  rotation-vibration  energy  levels  and 
wavefunctions  within  some  suitable  approximation.  This  approach  suffers 
from  the  fact  that  the  chosen  potential  usually  does  not  approximate 
the  true  potential  well  enough  to  predict  energy  levels  to  the  accuracy 
of  those  measured  by  microwave  and  high-resolution  infrared  techniques. 

Consequently,  most  workers  (refs.  1,2,15)  have  used  a third  approach 
in  which  the  vibration-rotation  Hamiltonian  is  a polynomial  expansion 
in  the  components  of  the  total  angular  momentum.  This  reduced  Hamiltonian 
has  the  advantage  that  it  does  not  require  any  information  about  the 
vibrational  motion  of  the  molecular  system,  as  only  the  angular  momentum 
operators  appear  in  it.  Instead,  it  contains  a large  number  of  rotation- 
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distortion  coefficients  that  must  be  extracted  from  the  observed  spectral 
line  data.  In  general,  an  asymmetric  rotor  such  as  water  requires 
(n  «■  1)  distortion  coefficients  to  each  order  n of  angular  momentum  oper- 
ator used  to  characterize  the  rotational  spectrum  within  a given  vibration- 
al band. 

Due  to  the  widespread  use  of  this  third  approach  one  would  natu- 
rally like  to  check  its  performance  within  some  model  system  where  analytic 
solutions  are  available.  Then,  by  comparing  the  predictions  of  the  ex- 
panded rotational  Hamiltonian  to  that  of  the  model,  a measure  of  its 
reliability  can  be  obtained.  In  particular,  the  following  five  features 
should  be  examined  in  detail. 

(1)  In  developing  a power  series  in  the  components  of  the 

total  angular  momentum,  what  physical  parameter  n is  this  expansion  based 
on,  i.e.,  is  n = rotational  energy/vibrational  energy? 

(2)  How  meaningful  is  this  characterization,  i.e.,  is  this  ap- 
proach limited  to  the  low-lying  J states  or  is  it  useful  for  all  rota- 
tional levels? 

(3)  Does  the  power  series  converge? 

(4)  Is  the  power  series  expansion  unique? 

(5)  To  what  extent  can  this  approach  predict  new  energy  levels? 

In  this  report,  these  questions  are  first  examined  within  the  context 
of  a rotating  diatomic  Morse  oscillator.  Then,  any  difficulties  that 
arise  in  this  simple  (diatomic)  system  can  be  expected  to  appear  to  a 
greater  degree  in  asymmetric  polyatomic  molecules  where  the  rotational 
structure  exhibits  a much  greater  degree  of  complexity  and  in  addition, 
possesses  a number  of  vibrational  degrees  of  freedom  as  well  as  coriolis 
forces . 

Also  note  that  for  highly  asymmetric  molecules  such  as  water, 
the  rotational  portion  of  the  molecular  wavefunction  depends  sensitively 
on  the  various  coefficients  that  appear  in  the  Watson  Hamiltonian.  As 
a consequence,  questions  (2)  through  (4)  are  of  critical  importance  for 
any  study  that  requires  accurate  wavefunctions . To  reinforce  the  results 
found  in  the  case  of  a rotating  Morse  oscillator,  some  of  the  same  tests 
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have  also  been  applied  to  the  water  molecule.  More  precisely,  standard 
techniques  were  used  to  construct  the  Watson  Hamiltonian  for  the  H2O 
molecule  and  then  questions  (2)  through  (S)  were  examined  for  this  case. 

Chapter  3 of  this  report  reviews  those  aspects  of  the  rotating 
Morse  oscillator  that  are  relevant  for  this  work,  discusses  the  physical 
parameter  that  arises  as  an  expansion  coefficient,  and  examines  the 
limitations  of  perturbation  expansion  techniques  as  they  appear  for  this 
system.  In  chapter  4,  the  predictions  of  the  reduced  Hamiltonian  are 
compared  to  the  analytic  results.  It  is  found  that  this  approach  is  not 
useful  in  discussing  the  higher  rotational  levels  (J  > 10)  of  light 
molecules  (4  'v  10  cm-1).  Chapter  5 is  devoted  to  the  numerical  analysis 
of  the  Watson  Hamiltonian  for  H2O  as  well  as  a discussion  of  the  findings. 
Finally,  chapter  6 contains  our  conclusions. 


SECTION  III 


THE  ROTATING  MORSE  OSCILLATOR 


Before  proceeding,  comments  are  presented  on  the  accuracy  of 
the  solutions  obtained  by  Pekeris  (ref.  16)  for  a rotating  Morse  os- 
cillator. Recall  that  Pekeris  neglected  all  terms  of  the  form 

AJ(J  * 1)0/  - l)3,  AJ{J  ♦ 1)0/  - 1)\... 

from  the  rotation- vibration  Hamiltonian.  Here  A is  the  rotational  con- 
stant, J the  rotational  quantum  number,  and 

y = exp[-a(r*  - rg)]  (1) 

In  equation  (1),  r is  the  equilibrium  bond  length,  r is  the  nuclear 
separation,  and  a characterizes  the  range  of  the  Morse  potential  t/^(r) , 
given  by 


£/(r)  = D(  1 - exp[-a(r  - r^)]}2  (2) 

Typically,  a = 1 X'1  and  D,  the  strength  of  the  potential,  is  about 
38,285  cm'1  for  H2  and  37,225  cm-1  for  HC1  (ref.  17).  To  estimate 
the  size  of  the  leading  term,  replace  the  quantity 


(r  - 


<ar 


(3) 


where  <ar2?  is  the  mean  square  deviation  of  the  nuclear  separation  from 
r • due  to  either  centrifugal  distortion  or  nuclear  vibration.  (As  a 
rough  estimate,  note  that  r^(n  = 1)  = 161/160  (r*  ) (n  = 0)  and 
^n-lJ-0~^n~0J-0  = ^00  Assume  that  the  stretching 

of  the  bond  from  its  equilibrium  length  varies  with  the  excitation 
energy  of  the  level  of  interest  (relative  to  the  ground  state),  then 


W|#V*“V 


5 x 10-2  ft 


and  taking  a = 1 X*1  find  that 


(4) 


AJ{J  * 1)Q,  - l)3  < 2 cm’1  (5) 

for  A = 10  ca*1.  By  way  of  comparison,  note  that  the  rotational 

energy  of  a Morse  oscillator  for  J = 40  is  about  13,000  cm'1  for  D = 
33,000  cm'1  and  for  A = 10  cm'1.  Thus,  we  may  safely  use  the  Pekeris 
solutions . 

Next,  the  energy  spectrum  of  a rotating  Morse  oscillator  is 
given  by  (ref.  16) 


AJ{J  + 1) 


1 


a2  = 


ar 


(ar  )2 
e 


(7c) 


If  one  expands  equation  (6)  in  powers  of  J (J  + 1),  the  first  two 
rotational  terms  are 


4-  ~ AJ(J  + 1)  - 

T1 


44  3 2 

J2(J  + 1) 


«0 


(8) 


where  ojq  = ap/2u)'i  is  the  harmonic  frequency,  and  u is  the  reduced 
mass.  For  the  expansion  to  be  physically  reasonable  one  would  conclude 
on  the  basis  of  equation  (8)  that 


(9) 


For  light  molecules  ojq  = 103  cm-1  and  A = 10  cm-1  so  that  J « 50, 
and  one  would  conclude  that  the  expansion  techniques  could  be  useful  for 
J < 15.  This,  however,  is  somewhat  optimistic.  An  examination  of 

equation  (6)  reveals  that  the  true  expansion  parameter  n is 


c2 


(10) 


For  ar  = 1,  find  that  a perturbation  expansion  is  valid  provided 


that 


J « 


D l"2 
2A 


(11) 
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For  D - 33  x 103  cm'1  and  A = 10  cm'1,  this  requires  that  J <<  40. 

Consequently,  perturbation  expansion  techniques  are  limited  to  J < 10. 

For  states  beyond  J = 10,  the  centripetal  distortion  energy  has  become 

so  large  that  a perturbation  expansion  based  on  the  rigid  rotor  approxi- 
mation requires  many  terms  to  converge.  As  a specific  example,  consider 
the  level  J = 20  for  which  n = 0.255.  To  ensure  an  accuracy  of  one 
part  in  106,  it  is  necessary  to  include  terms  of  order  n12  in  the  per- 


turbation expansion.  Now,  for  a diatomic  molecule,  the  rotational 
structure  is  characterized  by  only  one  angular  momentum  operator,  p2, 
and  such  an  analysis  is  possible.  On  the  other  hand,  the  situation  is 
virtually  hopeless  for  asymmetric  polyatomic  molecules.  There,  one  is 
faced  with  a variety  of  angular  momentum  operators  such  as  P2,  P2,  and 
p2  plus  many  combinations  such  as  ?4P2  and  so  on.  In  fact,  if  one  goes 
out  to  order  p2-n,  there  are  (2 n + 2)  distortion  coefficients.  As  a 
specific  example,  consider  water.  Note  that  water  has  three  different 
moments  of  inertia  of  which  the  two  largest  are  strongly  coupled  to  the 
bending  mode,  which  can  be  reasonably  well  approximated  as  a Morse 

oscillator  (ref.  5).  Then,  for  the  J = 16  level,  n - h,  and  this  requires 
terms  to  order  n20  in  order  to  characterize  all  of  the  energy  levels 
from  the  ground  state  to  J = 16,  to  one  part  in  106.  This,  in  turn, 
implies  that  one  requires  a perturbation  Hamiltonian,  which  includes 
terms  of  order  P40;  therefore,  there  are  several  hundred  unknown 
rotation-distortion  constants. 
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SECTION  IV 

NUMERICAL  COMPARISON 


This  section  examines  the  reliability  of  the  power  series 
expansion  approach  (as  generated  by  the  Watson  Hamiltonian  for  a dia- 
tomic molecule)  in  predicting  new  energy  levels  and,  in  addition,  tests 
its  uniqueness.  The  philosophy  here  is  to  regard  the  energy  level 
structure,  as  generated  by  equation  (6),  as  an  observed  spectrum.  Then, 
the  rotation-distortion  coefficients  that  appear  in  the  Watson  Hamiltonian 
are  fit  to  the  observed  low-lying  levels.  The  fit  should  be  to  the  low- 
lying  levels  to  ensure  that  the  moment  of  inertia  and  the  first  distortion 
coefficient  have  reasonable  values.  As  shall  be  seen,  if  one  fits  to 
the  higher  levels,  unphysical  values  are  obtained  for  the  distortion 
coefficients  and  rotation  constants. 

For  a diatomic  molecule,  the  Watson  Hamiltonian  reduces  to 

Jfp  = AXP2  + A2Pk  * A3P6  + ...  (12) 

where  P2  is  the  square  of  the  total  angular  momentum  and  the  A^  are  the 
rotation  distortion  coefficients  that  will  depend  on  the  vibrational 
band  but  should  not  vary  significantly  with  the  rotational  quantum  num- 
ber. Thus,  within  a given  vibrational  band,  the  rotational  energy  levels 
are 

Sj  = Eq  * AXJ(J  * 1)  + AzJ2{J  + l)2  + ...  (13) 

where  Eq  is  the  vibrational  energy. 

For  a specific  numerical  test,  the  following  values  have  been 
chosen  for  D,  o,  and  y: 

D = 3.3  x 104  cm'1;  a = 1.019  R"1 
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re  = 1-3264  %;  y = 1 amu. 

These  values  generate  a spectrum  that  roughly  describes  the  low-lying 
vibrational  states  of  the  bending  mode  of  H2O  (ref.  5). 

Recently,  De  Lucia,  Helminger,  Cook  and  Gordy  (ref.  1)  fit  15 
microwave  lines  of  H20  in  the  ground  vibrational  state  using  a pertur- 
bation expansion  that  included  terms  to  order  P10.  These  transitions 
involved  rotational  levels  as  low  as  J = 1 to  as  high  as  J = 10.  The 
equivalent  procedure  for  the  rotating  Morse  oscillator  would  be  to  fit 
those  levels  from  J = 1 to  J = 10  by  a perturbation  expansion  to  terms 
of  ?10.  Due  to  the  extreme  simplicity  of  the  rotational  structure  of 
diatomics,  this  means  that  the  following  spectrum 

Ej  =>  AXJ{J  ♦ 1)  ♦ A2Jz(J  + l)2  ♦ ...  + A5JS(J  + l)5  (14) 

should  be  fit  to  the  first  10  levels  generated  by  equation  (14).  (In 

equation  (14),  all  energies  are  measured  from  the  ground  rotational  state.) 
Using  a least  squares  fit,  the  values  in  table  1 were  obtained  for  the 
rotational-distortion  coefficients.  Finally,  note  that  the  various  co- 
efficients were  totally  uncorrelated;  consequently,  they  should  be 
highly  reliable.  (The  standard  deviations  for  the  rotation  coefficients 
B . and  C ■ (j  = 1-5)  are  zero  as  one  is  effectively  solving  for  five 

tJ  V 

unknowns  in  five  linear  equations.  The  standard  deviations  for  the  Aj 
and  D • will  be  nonzero.  However,  they  are  extremely  small  because 

V 

of  the  simple  character  of  the  energy  levels  and  the  fact  that  we  are 
fitting  to  only  10  levels.) 

A comparison  of  the  perturbation  expansion  (as  embodied  in 
equation  (14))  to  the  observed  energy  levels  is  exhibited  in  table  2 
for  the  first  30  rotational  levels.  An  examination  of  this  table  shows 
that  there  is  a significant  deviation  by  J - 13,  which  is  only  three 
states  and  656.20  cm"1  above  the  last  level  used  for  fitting  purposes. 

By  J = 20,  the  energy  level  structure  predicted  by  the  Watson  Hamiltonian 
deviated  by  more  than  325  cm"1  from  the  "observed"  value.  This  is  an 
8%  error  and  is  a signal  that  perturbation  expansion  techniques  are 
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breaking  down.  Another  indication  that  this  is  happening  can  be  seen 
by  comparing  the  first  term  in  equation  (14),  i.e.,  the  rigid  rotor 
energy,  to  the  second-order  term,  the  centripetal  distortion  energy. 

Thus , 

A\J(J  * 1)  - 4000  cm  * 

A2J2{_J  * l)2  * 500  cm"1 

for  J - 20.  When  this  occurs,  one  must  abandon  the  conceptual  picture 
of  a rigid  rotor  slightly  distorted  by  centripetal  forces.  Further 
evidence  that  the  rigid  rotor  picture  is  totally  inadequate  (and  there- 
fore the  perturbation  technique  is  inadequate)  is  dramatically  displayed 
by  the  rotational  levels  J = 22  and  J = 23,  where  the  predicted  energy 
levels  decrease  with  increasing  J.  Finally,  by  J = 27,  the  Watson 
Hamiltonian  predicts  negative  energy  values.  The  reason  for  this  dif- 
ficulty can  be  traced  to  the  fact  that  for  J = 27,  the  perturbation 
expansion  parameter  n = 0.45  and  many  more  terms  are  required  if  the 
series  is  to  converge.  Also  note  that  the  ratio  of  the  fifth  term  to 
the  first  is 

A^JSU  * l)5  « 2 

A\J(J  * 1) 

which  is  a clear-cut  sign  that  power  series  expansion  procedures  are 
essentially  useless. 

Table  1 

VALUES  FOR  ROTATIONAL  - DISTORTION  COEFFICIENTS 


Al 

= 

9.50271  x 

10°  cm'1 

a2 

3 

-1.45207 

X 

10* 3 cm"1 

A 3 

= 

-3.36340 

X 

10"7  cm  1 

Ak 

s 

+9. 11721 

X 

10  9 cm  1 

As 

= 

-4.31392 

X 

10  11  cm 

J 
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Table  2 

COMPARISON  OF  PERTURBATION  EXPANSION  TO  OBSERVED  ENERGY  LEVELS 


Observed 

Expansion 

Diff. 

(cm  J)* 

(cm'1) ‘ 

(cm  J) 

1 

19.036 

19.036 

0.0 

2 

57.073 

57.072 

-0.001 

3 

114.04 

114.039 

♦0.001 

4 

189.83 

189.832 

0.0 

S 

284.31 

284.312 

0.0 

6 

397.31 

397.306 

-0.005 

7 

528.61 

528.613 

0.0 

8 

678.00 

678.00 

0.0 

9 

845.20 

845.20 

0.0 

10 

1029.90 

1029.90 

0.0 

11 

1231.90 

1231.7 

-0.2 

12 

1450.7 

1450.03 

-0.67 

13 

1686.1 

1684.03 

-2.97 

14 

1937.6 

1932.31 

-5.29 

15 

2205.0 

2192.58 

-12.42 

16 

2487.6 

2661.11 

-17.49 

17 

2785.3 

2731.93 

-53.37 

18 

3097.5 

2995.68 

-101.82 

19 

3423.8 

3238.03 

-185.77 

20 

3763.8 

3437.54 

-326.26 

21 

4117.1 

3562.84 

22 

4483.2 

3568.85 

23 

4861.7 

3392.00 

24 

5252.2 

2944.12 

25 

5654.2 

2104.67 

26 

6067.4 

711.87 

27 

6491.4 

-1452.64 

28 

6925.7 

-4671.43 

29 

7370.0 

-9315.79 

30 

7823.9 

-15863.5 

*0 .01  cm  1 
tFor  p = 

'v  50  MHz  fitted  from  analytic 
0 band  only  expansion  to  F10 

solution 

i 


iL 


i, 


Ihe  rotation-distortion  coefficients  as  derived  from  different 
levels  are  not  unique.  Table  3 shows  a comparison  of  these  rotational- 
distortion  coefficients  used  to  fit  the  first  five  levels  (denoted 

bv  B ■)  versus  the  A..  Since  the  B-  were  fitted  to  the  first  five 

J J J 

rotational  levels,  they  are  nearly  equivalent  to  the  Taylor  series 
expansion  of  equation  (6)  to  P10  [i.e.,  to  e/5 (J  + l)5]  and,  consequently, 
they  are  the  exact  values  to  order  n6,  i.e.,  approximately  one  part 
in  1011.  Now,  Bi  = fi2/2 urg2  is  the  rotation  constant,  and  the 
corresponding  4-coefficient  deviates  from  the  correct  value  by 
three  parts  in  104.  This  is  far  more  serious  than  it  appears  because 
these  coefficients  are  often  quoted  to  parts  in  10®.  Furthermore,  this 
deviation  would  lead  to  an  incorrect  value  of  the  molecular  size  by 
parts  in  104.  Next,  consider  the  first  distortion  coefficient  where 
a deviation  of  over  2%  is  found.  Again,  this  is  serious  because  these 
coefficients  are  quoted  to  parts  in  104.  A further  examination  of  table 
3 reveals  that  there  is  almost  no  relationship  between  the  two  sets  of 
coefficients  beyond  j = 2.  Thus,  £5  is  almost  30  times  the  size  of 

A5,  and  they  are  of  opposite  signs.  One  might  argue  that  for  low- 
lying  levels  this  is  not  too  important.  Consider,  however,  J = 10. 

Then, 


A5J*(J  + l)5  = -0.430  cm" 


BSJS(J  + l)5  * +12.3  cm"1 


and  for  J - 20 


AsJ5iJ  + l)5 
B5J5(J  + l)5 


-430  cm  1 
+12,300  cm 


Note,  that  the  energy  associated  with  the  P10  term  is  about  three  times 
that  associated  with  the  rigid  rotor  energy.  The  implications  of  this 
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in  regard  to  calculating  asymmetric  rotor  wavefunctions  for  H2O  are 
significant  because  these  wavefunctions  depend  sensitively  on  the  rota- 
tion-distortion constant.  Also  note  that  the  large  discrepancy  between 
these  two  terms  implies  that  the  expansion  techniques  are  valid  only  for 
low-lying  states  where  the  lower  order  terms  dominate  the  level  structure. 
Finally,  in  table  4 are  listed  the  predictions  of  the  first  30  rota- 
tional levels  using  S1(  ....  S5. 

The  distortion  coefficients  are  very  sensitive  to  the  levels  that 
are  used  for  fitting  purposes.  Table  5 shows  the  distortion  coefficients 
that  arise  when  one  fits  the  J = 6 to  J = 10  levels.  These  coefficients 

are  denoted  by  C ..  Notice  that  significant  deviations  appear  immediately. 

C 

For  example,  the  rotation  constant  is  off  by  almost  one  part  in  104, 
and  is  larger,  not  smaller,  than  B^.  This  is  quite  serious  because  the 
first  coefficient  is  directly  related  to  the  moment  of  inertia.  Now, 
as  the  molecule  rotates  faster  and  faster,  it  gradually  opens  up,  i.e., 
it  stretches,  which  increases  its  moment  of  inertia.  Consequently,  if 
one  fits  the  rotation-distortion  coefficients  to  highly  excited  rotational 
states,  a larger  moment  of  inertia  would  be  expected  and,  therefore,  a 
smaller  rotation  constant.  Instead,  we  find  that 


The  first  distortion  coefficient  is  off  by  2%.  Also,  the  last  two  co- 
efficients are  between  one  and  two  orders  of  magnitude  too  small.  Thus, 
one  may  draw  the  conclusion  that  the  distortion  coefficients  are  not  only 
not  unique,  but  they  vary  in  a sensitive  and  nonphysical  manner  with 
whichever  states  are  used  for  fitting  purposes.  Also  note  that  the 
large  difference  in  size  of  the  higher  order  distortion  constants  casts 
doubt  on  their  physical  meaning.  Finally,  the  resultant  spectrum  is 
displayed  in  table  6. 
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Table  3 

COMPARISON  OF  ROTATIONAL  - DISTORTION  COEFFICIENTS  VS.  A ■ 

J 

J B 


1 

9.52101 

9.52071 

+0.0030 

3.14  x 10 

2 

-1.48844 

X 

10*3 

-1.45207 

X 

10'3 

-0.0313  x 10'3 

2.1  x 10' 

3 

+1 . 58233 

X 

10'6 

-3.36340 

X 

10'7 

+1.94573  x 10'6 

5.7 

4 

-6.94441 

X 

10'8 

9.11721 

X 

10'9 

-7.856  x 10"8 

8.6 

5 

1.24007 

X 

10'9 

-4.31392 

X 

10"u 

+1.2833  x 10*9 

29 

Table  4 

PREDICTIONS  OF  FIRST  30  ROTATIONAL  LEVELS 


J_  Expansion  (cm  *)  Observed  (cm 


1 

1.90361 

X 

101 

19.036 

2 

5.70729 

X 

101 

57.073 

3 

1.14040 

X 

102 

114.04 

4 

1.89832 

X 

102 

189.83 

5 

2.84211 

X 

102 

284.31 

6 

3.97327 

X 

102 

397.31 

7 

5.28799 

X 

102 

528.61 

8 

6.78941 

X 

102 

678.00 

9 

8.48787 

X 

102 

845.20 

10 

1.04126 

X 

103 

1029.90 

11 

1.26316 

X 

103 

1231.90 

12 

1.52860 

X 

103 

1450.7 

13 

1.86463 

X 

103 

1686.1 

14 

2.32000 

X 

103 

1937.6 

15 

2.97844 

X 

103 

2205.0 

16 

3.97787 

X 

103 

2487.6 

17 

5.53791 

X 

103 

2785.3 

18 

7.99781 

X 

103 

3097.5 

19 

1.18682 

X 

104 

3423.8 

20 

1.78999 

X 

104 

3763.8 

21 

2.71751 

X 

104 

4117.1 

22 

4.12240 

X 

104 

4483.2 

23 

6.21757 

X 

104 

4861.7 

24 

9.29479 

X 

104 

5252.2 

25 

1.37484 

X 

105 

5654.2 

26 

2.01048 

X 

105 

6067.4 

27 

2.90584 

X 

105 

6491.4 

28 

4.15160 

X 

10s 

6925.7 

29 

5.86494 

X 

10s 

7370.0 

30 

8.19598 

X 

105 

7823.9 

B constants  are  : . 9.52101,  -1.48444  x 10~3,  1.58233  x 10’6, 
-6.94441  x 10‘6,  and  1.24007  x 10‘9 
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Table  5 

DISTORTION  LEVELS  FROM  FITTING  J = 6 TO  J = 10  LEVELS 


J 

1 

2 

3 

4 

5 


B .(cm'1) 

c 

9.52101 

-1.4844  x 10-3 
1.5823  x 10*6 
-6.9444  x 10-8 
1.2401  x 10-9 


(cm-1) 

9.52187 

-1.52611  x 10"3 
+1.40507  x 10’6 
-9.0708  x IQ'9 


-11 


2.83792  x 10 


Table  6 


RESULTANT  SPECTRUM 


<1  Expansion(cm~*)  Observed  (cm-1) 


1 

1.90375 

X 

101 

19.036 

2 

5.70768 

X 

101 

57.073 

3 

1.14047 

X 

102 

114.04 

4 

1.89843 

X 

102 

189.83 

5 

2.84324 

X 

102 

284.31 

6 

3.97320 

X 

102 

397.31 

7 

5.28626 

X 

102 

528.61 

8 

6.78011 

X 

102 

678.00 

9 

8.45213 

X 

102 

845.20 

10 

1.02995 

X 

103 

1029.90 

11 

1.23190 

X 

103 

1231.90 

12 

1.45074 

X 

103 

1450.7 

13 

1.68611 

X 

103 

1686.1 

14 

1.93765 

X 

103 

1937.6 

15 

2.20496 

X 

103 

2205.0 

16 

2.48764 

X 

103 

2487.6 

17 

2.78529 

X 

io? 

2785.3 

18 

3.09748 

X 

iO3 

3097.5 

19 

2.42379 

X 

103 

3423.8 

20 

3.76380 

X 

IO3 

3763.8 

21 

4.11709 

X 

IO3 

4117.1 

22 

4.48328 

X 

10  3 
10, 

4483.2 

23 

4.86198 

X 

4861.7 

24 

5.25286 

X 

io3 

5252.2 

25 

5.65566 

X 

103 

5654.2 

26 

6.07018 

X 

103 

6067.4 

27 

6.49635 

X 

10 

10 

10? 

6491.4 

28 

6.93423 

X 

6925.7 

29 

7.38407 

X 

7370.0 

30 

2.84638 

X 

IO3 

7823.9 

31 

8.32198 

X 

10  3 

8287.1 

32 

8.81208 

X 

10 

10? 

8759.2 

33 

9.31837 

X 

9239.9 

34 

9.84318 

X 

10 

9728.9 

35 

1.03896 

X 

IO* 

10226 

36 

1.09615 

X 

io4 

10730 

37 

1.15641 

X 

104 

11242 

38 

1.22038 

X 

ioi 

11762 

39 

1.28886 

X 

io4 

12288 

40 

1.36284 

X 

io4 

1.2821  x IO4 

1 


j. 


\ 
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To  show  how  misleading  perturbation  expansion  techniques  can  be, 
table  7 shows  the  results  of  fitting  the  J = 11  to  J = 20 
rotational  levels  in  the  ground  vibrational  state.  The  rotational  con- 
stants, denoted  by  D.,  were  found  by  a least  squares  fit  displayed 
in  table  8. 

At  first  glance  it  would  appear  that  this  fit  is  far  better  than 
the  previous  one  in  which  the  spectrum  was  fit  to  the  first  10  rotational 
levels.  This,  however,  is  not  the  case.  An  examination  of  table  7 
reveals  that 


which  is  opposite  to  what  would  be  expected.  Another  disturbing  feature 
of  table  8 is  that  the  last  three  coefficients  are  more  than  an  order 
of  magnitude  smaller  than  the  B values,  and  the  last  coefficient  is 
three  orders  of  magnitude  smaller.  These  features  reinforce  the  con- 
clusion that  the  higher-order  rotation  distortion  coefficients,  generated 
by  the  Watson  Hamiltonian,  are  devoid  of  all  physical  significance  for 
light  molecular  systems  that  are  excited  to  high  rotational  states.  They 
are,  in  fact,  just  a sequence  of  numbers  that  are  chosen  to  give  a fit 
of  a power  seri'es  expansion  of  the  form 


E = l CnJn(J  * l)n 
n 

to  the  rotational  levels.  It  does  this  by  vastly  decreasing  the  magni- 
tude of  the  higher-order  rotation-distortion  coefficients,  as  these  co- 
efficients have  the  greatest  weight  associated  with  them  for  high  J. 
Furthermore,  note  that  if  one  examines  even  higher  lying  rotational 
states,  deviations  again  arise,  e.g.,  by  J - 40,  the  discrepancy  is  over 
800  cm  * . Also,  note  that  the  low-lying  levels  are  in  reasonable 
agreement.  This  arises  because  the  energy  of  these  states  is  determined 
by  the  first  few  rotation-distortion  coefficients,  which  are  always 
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correct  within  parts  in  104.  Finally,  note  that  the  J = 8,9,10 
states  are  in  agreement  because  they  lie  near  the  levels  used  for  the 
purpose  of  fitting. 

In  summary,  within  the  context  of  the  rotating 
Morse  oscillator,  the  perturbation  expansion  technique,  when  applied  to 
highly  excited  rotational  states  of  light  molecules,  is  devoid  of  any 
physical  significance.  This  arises  from  the  fact  that  the  perturbation 
expansion  parameter  n = 2AJ(J  + 1 )/D  becomes  large  (-**  to  h)  for 
light  molecules  near  J = 10.  As  a consequence,  the  series  is  a slowly 

convergent  one  requiring  many  terms.  Attempts  to  fit  the  energy 
spectrum  with  just  a few  terms  can  be  done  only  by  significantly 
changing  the  distortion  coefficients  so  that  they  bear  little  resem- 
blance to  reality.  This  is  borne  out  by  the  fact  that  if  one  examines 
even  higher  levels,  the  predicted  spectrum  is  still  incorrect  by  hundreds 
of  wave  numbers.  Thus,  the  perturbation  expansion  technique  cannot 
predict  new  energy  levels.  Furthermore,  it  is  clear  that  one  can  charac- 
terize only  the  low-lying  levels  (where  the  rotational  distortion  energy 
is  small  compared  to  the  total  rotational  energy)  by  means  of  the  Watson 
Hamiltonian  approach.  Also  note  that  the  results  are  not  unique,  i.e., 
if  one  fits  the  J = 1 to  J = 5 levels  or  the  J = 6 to  10,  different 
distortion  coefficients  are  found. 
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Table  7 


RESULTS  OF  FITTING  THE 

J 

= 11 

TO  / = 

20  ROTATIONAL  LEVELS 

J_ 

Expansion 

Cent”1) 

Observed  (cm”1; 

1 

1.90376 

X 

101 

19.036 

2 

5.70766 

X 

10 

57.073 

3 

1.14045 

X 

102 

114.04 

4 

1.89837 

X 

102 

189.83 

5 

2.84314 

X 

1°2 

284.31 

6 

3.97306 

X 

10 

397.31 

7 

5.28612 

X 

102 

528.61 

8 

6.77999 

X 

102 

678.00 

9 

8.45204 

X 

10 

10 

845.20 

10 

1.02994 

X 

1029.90 

11 

1.23191 

X 

103 

103 

1231.90 

12 

1.45086 

X 

1450.7 

13 

1.68662 

X 

10  3 

1686.1 

14 

1.93927 

X 

10 

1937.6 

15 

2.20929 

X 

10 

2205.0 

16 

2.49795 

X 

10 

2487.6 

17 

2.80771 

X 

10 

2785.3 

18 

3.14295 

X 

10  3 
103 

3097.5 

19 

3.51088 

X 

3423.8 

20 

3.92288 

X 

103 

2763.8 

21 

4.39623 

X 

103 

4117.1 

22 

4.95641 

X 

10 

4483.2 

23 

5.64011 

X 

10 

4861.7 

24 

6.49903 

X 

10 

5252.2 

25 

7.60476 

X 

10 

5654.2 

26 

9.05482 

X 

104 

6067.4 

27 

1.09801 

X 

io4 

6491.4 

28 

1.35543 

X 

104 

6925.7 

29 

1.70049 

X 

104 

7370.0 

30 

2.16265 

X 

104 

7823.9 

31 

2.77978 

X 

104 

8287.1 

32 

3.60001 

X 

104 

8759.2 

33 

4.68405 

X 

i°44 

9239.9 

34 

6.10790 

X 

10 

9728.9 

35 

7.96599 

X 

105 

10226 

36 

1.03749 

X 

10 

10730 

37 

1.34774 

X 

io! 

11242 

38 

1.74478 

X 

10c 

11762 

39 

2.24972 

X 

5 

10c 

12288 

40 

2.88803 

X 

10s 

1.2821  x H 

L 
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Table  8 

FOR  ROTATIONAL  CONSTANTS 

D .(cm' 

J 

9.52170  -0.00069 

-1.4871  x 10“3  +0.003  x 10'3 

+4.4446  x 10'7  1.1378  x 10'6 

-2.5619  x 10'9  -6.6832  x 10'9 

+1.6066  x 10-11  1.2384  x 10'10 


(B.  - DOC  cm'1) 
sZ sZ_ 


SECTION  V 


NUMERICAL  ANALYSIS  FOR  THE  WATER  MOLECULE 

In  the  previous  two  sections  the  validity  of  power  series  ex- 
pansion techniques  was  examined  for  a model  system.  In  this  section, 
similar  tests  are  applied  to  a real  polyatomic  asymmetric  molecule, 
namely  H2O.  Owing  to  the  mathematical  complexity  of  extreme  asymmetric 
rotors,  no  analytic  work  appears  possible.  Hence,  the  convergence, 
uniqueness,  and  predictability  of  power  series  expansion  techniques 
have  only  been  examined  numerically. 

There  are  several  techniques  used  in  analyzing  the  Watson 
rotation  Hamiltonian.  The  simplest  consists  of  a least  squares  analysis 
to  derive  the  various  rotation-distortion  coefficients.  The  resulting 
Hamiltonian  is  then  used  to  derive  energy  levels  within  a rigid  asym- 
metric rotor  basis  (as  determined  by  the  A,  B,  and  C rotation  constants), 
and  no  attempt  is  made  to  diagonalize  the  rotation  Hamiltonian.  A 
second  approach,  perhaps  the  most  widely  used,  involves  a matrix  diag- 
onalization  in  addition  to  a least  squares  analysis.  In  the  appendix, 
it  is  shown  that  a matrix  diagonali zation  does  not  give  rise  to  signifi- 
cant changes  from  a linearized  analysis.  Finally,  a third  technique, 
used  by  De  Lucia,  et  al.,  is  designed  to  extract  out  higher  order  dis- 
tortion energies  from  the  distortion  coefficients.  This  approach  is 
discussed  in  reference  1 and  we  have  found  it  to  be  far  superior  to  the 
first  two.  Nevertheless,  it  has  a number  of  defects  insofar  as  H20  is 
concerned.  In  general,  the  numerical  results  for  H20  confirm  the  findings 
for  the  rotating  Morse  oscillator;  namely,  one  cannot  deduce  an  adequate 
physical  picture  of  the  molecule  without  explicitly  including  the  coupling 
of  the  vibrational  motion  of  the  nuclei  to  the  molecular  rotation. 

CONVERGENCE  OF  THE  POWER  SERIES 

As  a first  check  on  the  usefulness  and  validity  of  power  series 
expansions,  we  have  examined  the  rate  at  which  these  series  converge. 


This  is  an  essential  test  because  if  the  perturbation  expansion  (about 
a rigid  asymmetric  rotor)  is  to  be  meaningful  it  must  converge  rapidly. 
Although  any  criteria  devised  for  convergence  are  to  some  degree 
arbitrary,  each  term  in  the  series  should,  in  general,  be  about  one 
order  of  magnitude  smaller  than  the  term  proceeding  it.  Then,  for  a 
rotational  level  that  lies  about  one  thousand  wave  numbers  above  the 
rotational  ground  state,  the  sum  of  the  energies  of  all  the  tenth-order 
terms  (to  date,  these  are  the  highest  order  terms  used)  should  be  on  the 
order  of  one -tenth  of  a wave  number.  This,  in  turn,  ensures  an  accuracy 
of  between  one-hundredth  and  one-tenth  of  a wave  number. 

As  a specific  example,  the  rotational  Hamiltonian  given  by 
De  Lucia,  et  al.,  (refs.  1,2)  was  used  to  calculate  the  energy  contribution 
that  arises  in  each  order.  The  various  rotation  and  rotation-distortion 
constants  that  appear  were  obtained  from  microwave  and  infrared  data 
involving  the  low-lying  rotational  states.  Consequently,  it  should  be 
the  most  accurate  Hamiltonian  available,  and  any  convergence  problems 
that  arise  here  should  also  appear  with  other  Watson  Hamiltonians. 

One  might  argue  that  these  constants  should  be  extracted  from  data  that 
include  higher  rotational  levels.  However,  the  resulting  Hamiltonian 
is  less  accurate,  and  this  point  is  more  thoroughly  discussed  in  the 
next  subsection. 

The  convergence  calculation  was  carried  out  in  a rigid  asymmetric 
rotor  basis.  More  precisely,  we  evaluated  the  diagonal  contribution 
of  each  term  of  the  Hamiltonian  of  De  Lucia,  et  al.,  using  the  eigen- 
functions of  an  asymmetric  rotor  with  rotation  constants 

A = 27.8806  cm-1 

B = 14.5217  cm'1 

C = 9.27769cm'1. 


Note:  References  1 and  2 are  concerned  only  with  the  low  energy  states. 

The  results  of  reference  12  are  used  only  to  illustrate  the  point  that 
many  higher  order  terms  (P12,  P14,  etc.)  must  be  included  in  the  Watson 
Hamiltonian  to  adequately  describe  the  water  molecule  for  the  higher 
order  state. 
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This  corresponds  to  a linear  analysis  instead  of  a matrix  diagonalization. 
However,  the  difference  between  the  energies  obtained  using  a linear 
analysis  and  matrix  diagonalized  energies  is  generally  less  than  0.05  cm"1 
for  all  states  through  J equals  10  and  under  five  wave  numbers  for  all 
rotational  states  through  J equals  16. 

Table  9 lists  all  the  various  contributions  to  the 
energy  for  all  rotational  levels  through  J equals  10.  Here,  En  refers 
to  the  rigid  rotor  energy,  Ek,  ff6,  E 9,  and  Ey  0 to  the  fourth,  sixth, 
eighth,  and  tenth  order  perturbation  energies,  and  E^  is  the  total  cal- 
culated energy.  The  next  to  last  column  refers  to  the  experimental 
value  of  the  energy  as  deduced  from  spectroscopic  data  (ref.  18),  and  A E is 
equal  to  the  difference  between  the  matrix  diagonalized  energy  and  the 
linearized  energy.  All  energies  are  in  units  of  wave  numbers. 

Recall  that  our  analysis  in  chapter  2 indicates  that 
the  perturbation  expansion  parameter  n was 

n = 2AJ{J  + 1 )/D  (5.1) 

where  AJ(J  * 1)  is  the  rigid  rotor  energy  and  D is  the  strength  of  the 
Morse  potential.  As  mentioned  in  chapter  2,  the  bending  mode  of  H2O 
can  be  approximated  by  a Morse  oscillator,  and  it  is  reasonable  to 
adopt  equation  (5.1)  as  the  expansion  parameter  for  this  case,  although  the 
coefficient  (i.e.,  the  factor  of  2)  is  undoubtedly  different.  Con- 
sequently, take 


nH20 


£.  (000) ground  state  asymmetric  rotor  energy 

33,000  cm-1 


(5.2) 


as  the  expansion  parameter.  Here,  f is  a numerical  factor  that  lies 
between  1 and  5,  say,  and  we  have  taken  D equal  to  3.3  x 104  cm"1. 

Then,  for  C - 3 and  a rigid  rotor  energy  of  2750  cm'1,  the  fifth-order 
energy  would  be  approximately  a few  wave  numbers.  For  H20,  convergence 
difficulties  would  arise  at  the  IO.7  level. 

An  examination  of  table  9 reveals  that  the  tenth  oTder  pertur- 
bation energies  of  all  the  rotational  levels  through  J = 6 are  one- 
tenth  of  a wave  number  or  less.  Beyond  J = 6,  there  are  a number  of 


r 


levels  that  do  not  converge  as  rapidly.  For  example,  consider  the  I7 
and  76  states,  where  the  terms  decrease  by  a factor  of  5 only,  and 
the  last  term  is  0.5  cm’1.  This  is  approximately  2 orders  of  magnitude 
greater  than  A E and,  consequently,  cannot  be  attributed  to  the  fact  that 

only  a linearized  calculation  was  carried  out.  Instead,  it  implies  that 
at  least  for  this  level,  twelveth-order  terms  should  be  included. 

Deviations  for  such  low  levels  might  be  attributed  to  the 

fact  that  the  microwave  data  used  by  De  Lucia,  et  al.,  contained  information 
on  only  the  low  t states  and,  consequently,  difficulties  might  be  expected 
for  high  t states.  Another  possibility  is  that  the  effective  expansion 
parameter  given  by  equation  (5.2)  is  sufficiently  large  to  prevent  convergence. 
This  would  be  the  case  for  C = 4.  The  situation  for  J - 8 is  similar, 

although  somewhat  worse.  Here,  the  85  and  85  states  have  tenth-order 
energies  that  are  approximately  half  a wave  number,  and  the  87  and 
80  levels  are  two  wave  numbers.  In  both  cases,  the  convergence  is  slow, 
which  implies  that  P12  and,  possibly,  P14  terms  should  be  included  in  the 
rotational  Hamiltonian. 

Further  examination  of  table  9 reveals  that  convergence  problems 
arise  for  the  94  through  the  99  states.  Also,  the  IO5  through  10*  0 
states  all  have  tenth-order  energies  that  are  greater  than  one  wave  number. 

An  adequate  description  of  these  levels  would  probably  require  P12  and 
P14  order  terms  and  in  some  cases,  possibly,  P16  due  to  the  slow  convergence. 

One  might  argue  that  since  these  levels  differ  by  only  a few 
tenths  of  a wave  number  from  the  observed  energies,  a tenth-order  Hamil- 
tonian is  actually  sufficient.  However,  this  argument  is  misleading. 

The  higher-order  energies  of  those  states  that  are  well  described  by  this 
Hamiltonian  (i.e.,  all  levels  through  J = 6 and  all  the  small  t levels) 

have  tenth-order  energies  that  are  vanishingly  small.  Also  note  that 
these  are  the  levels  that  are  involved  in  the  microwave  data,  which  con- 
tain no  information  on  the  high  J T states.  The  total  energies  of  the 
states  agree  reasonably  well  with  experiment  because  these  levels  lie 
only  a few  hundred  wave  numbers  away  from  the  measured  ones. 

Table  10  lists  levels  that  lie  between  J = 11  and 

J = 16  in  which  the  fifth-order  term  is  greater  than  1 cm-1.  In  a 

few  isolated  cases,  this  difficulty  might  be  attributed  to  the  fact  that 
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a linear  analysis  was  performed  instead  of  an  exact  matrix  diagonali- 
zation.  However,  in  the  vast  majority  of  cases,  this  is  not  true  as  is 
borne  out  by  the  analysis  in  the  appendix.  In  general,  the  matrix 
diagonalization  gives  rise  to  a correction  of  under  1 cm"1. 

Consider,  first,  the  J = il  levels.  For  the  high  1 states,  i.e., 
r = 7 to  11,  the  convergence  factor  for  the  last  few  terms  is  about  3. 

Since  Pio  = 45  cm"1  for  r = 11,10  and  Ejq  = 18  cm'1  for  r = 8,9, 

one  would  have  to  include  terms  of  order  P18  to  describe  all  the  J = 11 

levels  and  retain  spectroscopic  accuracy  for  the  microwave  spectrum. 

Another  interesting  feature  displayed  in  table  10  is  that  although  the  total 
energy  of  the  various  states  is  often  incorrect  by  only  a few  tenths  of 
a wave  number  from  the  observed  energy,  the  series  do  not  con- 
verge. Presumably,  the  higher  order  terms  alternate  in  sign  so  that 
their  effect  on  E ^ is  small  or  higher  order  effects  are  being  included  in 
the  lower  coefficients.  For  example,  the  t = 6 and  7 levels  could 

require  an  E\2  as  well  as  an  Pi 4 although  | E^  - P0BS|  = 0.26  cm-1. 

(Although  the  convergence  is  poor  for  the  high  t states  of  the  J = 11 

levels,  the  physical  picture  of  a rigid  rotor  slightly  perturbed  by 
distortion  is  still  viable.) 

For  J = 12  to  14  the  situation  is  similar  except  that  the  con- 

vergence is  now  about  2 for  the  high  x states,  and  the  Pjq  terms  have 
grown  in  both  absolute  magnitude  as  well  as  relative  size  in  comparison 
to  E^.  For  these  states,  the  physical  picture  of  an  asymmetric  rotor 
slightly  perturbed  by  the  rotation-vibration  interaction  is  no  longer 
valid.  This  is  reflected  by  the  fact  that  for  J ' = 13i3,  Ejo  is  5% 

of  the  rigid  rotation  energy.  Owing  to  the  slow  rate  of  convergence  of 
the  rigid-rotor  model,  it  is  virtually  impossible  to  describe  these 
states  accurately  and  still  retain  a detailed  picture  of  the  micro- 
wave  spectrum. 

For  the  high-lying  t levels  of  the  states  with  J = 15,  there 

is  no  convergence.  In  some  situations,  higher  order  rotation-distortion 
terms  are  larger  in  magnitude  than  some  of  the  lower  order  ones. 

Another  feature  is  that  states  with  x < 0 are  not  converging  rapidly. 

Such  states  do  not  have  a great  deal  of  rotation-vibration  interaction 
and  probably  could  be  correctly  described  if  P12  and  P14  terms  were 
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included.  Furthermore,  for  J = 1(\  the  situation  is  essentially  the 

same.  However,  a new  feature  has  appeared.  States  that  lie  lower  in 
energy  than  the  J = 15  states  no  longer  satisfy  the  criteria.  Thus, 

the  t = -10  state  lies  lower  in  energy  than  the  15 _ 3 level,  and  it 

does  not  meet  our  convergence  requirements.  Finally,  on  the  basis  of 
the  trend  in  able  10,  it  appears  likely  that  for  J > 17,  convergence 

difficulties  will  arise  for  all  x-states. 

UNIQUENESS 

Next,  we  examine  the  dependence  (if  any)  of  the  rotation  and 
higher  order  distortion  constants  on  the  particular  set  of  energy  levels 
used  to  obtain  them.  In  this  test,  no  attempt  was  made  to  remove  the 
effects  of  the  rotation-vibration  interaction  from  these  coefficients 
(as  was  done  in  reference  12) . Instead,  a least  squares  analysis  was  carried 
out  on  the  rotation  and  higher  order  distortion  constants  of  the  Hamiltonian 


Jf  = * C)  P2  + [A  - H (fl  + C)](P2-b2P 2) 

z p - 

- V‘  - V/2pI  - V - 2{jp2pi  -Vi2?2  * ?2f2> 

♦ V6  * Vhpi  * vAI  * VI  * 2V4pI 

+ + p2p2)  + ^(PS4P2  + P2P4) 

* £VY  * W2pI  * VI  * vfpl  * p5pI> 

+ P^P10  + PX(P8P2  + P2P8), 


(5.3) 


to  ensure  that  the  energy  levels  of  a selected  set  of  states  of  equation 
(5.3)  agree  with  available  experimental  data.  This  approach  has  been  used 
by  a number  of  workers,  and  results  of  this  work  reveal  some  of  the 
limitations  inherent  in  it. 


L 


29 


CONVERGENCE  01-  THE  POWER  SERIES 
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ENERGY  OF 

VARIOUS 

TERMS  OF 

THE  HAMILTONIAN 

J 

T 

er 

*6 

ElO 

et 

Observed 

11 

11 

3506.95 

-404.92 

186.86 

-114.63 

45.46 

3219.71 

3216.6 

11 

10 

3506.95 

-404.92 

186.86 

-114.63 

45.46 

3219.71 

3216.6 

11 

9 

317.43 

-270.26 

102.69 

-53.19 

17.99 

2973.65 

2973.07 

11 

3 

3176.43 

-270.26 

102.69 

-53.19 

17.99 

2973.65 

2973.07 

11 

7 

2878.15 

-173.62 

52.22 

-22.41 

6.13 

2740.47 

2740.73 

11 

6 

2878. 15 

-173.62 

52.22 

-22.41 

6.13 

2740.47 

2740.73 

11 

5 

2612.45 

-107.53 

23.89 

-8.15 

1.59 

2522.25 

2522.46 

11 

4 

2612.45 

-107.53 

23.89 

-8.15 

1.59 

2522.25 

2522.46 

12 

12 

4160.32 

-573.97 

315.26 

-230.13 

108.43 

3779.92 

3767.1 

12 

11 

4160.32 

-573.97 

315.26 

-230.13 

108.43 

3779.92 

3767.1 

12 

10 

3798.26 

-396.77 

182.62 

-114.07 

46.51 

3516.56 

3312.8 

12 

9 

3798.26 

-396.77 

182.62 

-114.07 

46.51 

3516.56 

3512.8 

12 

8 

3468.41 

-265.59 

99.29 

-52.25 

17.72 

3267.56 

3267.2 

12 

7 

3468.41 

-265.59 

99.29 

-52.25 

17.72 

3267.56 

3267.2 

12 

6 

3171.04 

-172.18 

49.83 

-21.38 

5.55 

3032.65 

3033.17 

12 

5 

3171.04 

-172.18 

49.63 

-21.38 

5.55 

3032.65 

3033.17 

12 

4 

2906.64 

-109.09 

22.05 

-7.21 

1.09 

2813.48 

2813.94 

12 

3 

2906.64 

-109.09 

22.05 

-7.21 

1.09 

2813.48 

2813.94 

13 

13 

4869.44 

-791.13 

510.03 

-436.87 

241.28 

4392.76 

_ _ 

13 

12 

4869.44 

-791.13 

510.03 

-436.87 

241.28 

4392.76 

-- 

13 

11 

4475.86 

-563.24 

308.78 

-228.93 

107.55 

4103.21 

-- 

13 

10 

4465.86 

-563.24 

308.78 

-228.93 

107.55 

4103.21 

-- 

15 

9 

4114.44 

- 90.17 

177.31 

-112.33 

46.13 

3835.38 

-- 

13 

8 

4114.44 

-390.17 

177.31 

-112.33 

46.14 

3835.38 

• - . 

13 

7 

3785.43 

-262.84 

95.11 

-50.40 

16.57 

3584.0 

3583.87 

13 

6 

3785.43 

-262.84 

95.11 

-50.40 

16.57 

3584.0 

3583.87 

13 

5 

3489.24 

-173.02 

46.52 

-19.68 

4.47 

3347.52 

3348.2 

13 

4 

3489.24 

-175.02 

46.52 

-19.68 

4.47 

3347.52 

3348.2 

14 

14 

5645.45 

-1064.77 

796.17 

-790.83 

506.0 

5080.91 

_ _ 

14 

13 

5634.34 

-1064.77 

796.17 

-790.83 

506.0 

5080.91 

14 

12 

5209.21 

-777.36 

500.45 

-434.58 

246.08 

4743'.  81 

14 

11 

5209.21 

-777. 36 

500.45 

-434.58 

246.08 

4743.81 



14 

10 

4816.24 

-554.30 

300.80 

-22S.88 

110.36 

4447.22 

14 

9 

4816.24 

-554.30 

300.80 

-225.88 

110.36 

4447.22 

14 

8 

4455.63 

-385.24 

170.88 

-109.13 

44.04 

4175.67 

_ - 

14 

7 

4455.63 

-385.24 

170.88 

-109.13 

44.04 

4175.67 

-- 

14 

6 

4127.09 

-262.64 

90.14 

-47.42 

14.40 

3922.17 

— 

14 

S 

4127.69 

-262.64 

90.14 

-47.42 

14.40 

3922.17 

_ - 

14 

4 

3832.98 

-176.80 

42.95 

-17.14 

2.80 

3684.78 

-- 

14 

3 

3832.97 

-176.80 

42.95 

-17.14 

2.80 

3684.47 

— 

34 
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Table  10  (Continued) 


J 

T 

** 

*6 

*8 

*10 

Observed 

15 

15 

6454.99 

-1403.93 

1205.16 

-3174.06 

1008.29 

5890.05 

15 

14 

6454.99 

-1403.93 

1205.16 

-1324.06 

1008.29 

5890.05 

-- 

15 

13 

5998.33 

-1047.42 

782.41 

-786.72 

515.43 

5462.03 

-- 

15 

12 

5998.33 

-1047.42 

782.41 

-786.72 

515.43 

5462.03 

-- 

15 

11 

5573.81 

-765.58 

488.83 

-429.43 

245.97 

5113.59 

15 

10 

5573.58 

-765.58 

488.83 

-429.43 

245.97 

5113.59 

-- 

15 

9 

5181.60 

-547.77 

291.24 

-220.53 

106.75 

4811.28 

-- 

15 

8 

5181.60 

-547.77 

291.24 

-220.53 

106.75 

4811.28 

-- 

15 

7 

4821.98 

-384.14 

163.28 

-104.13 

39.93 

4536.93 

-- 

15 

6 

4821.98 

-384. 14 

163.28 

-104.13 

39.93 

4536.93 

-- 

15 

5 

4495.38 

-265.69 

84.43 

-43.08 

11.02 

4282.06 

-- 

IS 

4 

4495.38 

-265.69 

84.43 

-43.08 

11.02 

4282.06 

-- 

15 

3 

4202.53 

-184.29 

39.03 

-13.61 

4.31 

4044.10 

4045.8 

15 

2 

4202.50 

-184.28 

39.03 

-13.61 

4.30 

4044.06 

4045.8 

15 

1 

3945.04 

-132.88 

15.57 

-1.13 

-2.25 

3824.35 

2836.1 

15 

0 

3944.38 

-132.65 

15.48 

-1.08 

-2.27 

3823.86 

3824.8 

15 

-1 

3728.41 

-106.86 

6.24 

2.75 

-2.05 

3628.5 

3628.7 

15 

-2 

3721.27 

-103.65 

5.23 

3.10 

-2.20 

2623.75 

3624.0 

15 

-3 

3568.22 

-103.60 

5.88 

2.69 

-1.19 

3472.0 

3473.0 

16 

16 

7331.41 

-1818.30 

1775.99 

-2303.69 

1921.77 

6907.18 



16 

15 

7331.41 

-1818.30 

1775.99 

-2303.69 

1921.77 

6907.18 

-- 

16 

14 

6843.21 

-382.43 

1185.89 

-367.08 

1025.98 

6305.57 

-- 

16 

15 

6843.21 

-1382.43 

1185.89 

-1367.08 

1025.98 

6305.37 

-- 

16 

12 

6387.14 

-1032.31 

765.93 

-778.33 

516.30 

5858.73 

-- 

16 

11 

6387.14 

-1032.31 

765.93 

-778.33 

516.30 

5858.73 

-- 

16 

10 

5963.35 

-756.50 

475.03 

-420.79 

240.09 

5501.19 

-- 

16 

9 

5963.35 

-756.50 

475.03 

-420.79 

240.09 

5501.19 

-- 

16 

8 

5572.08 

-544.36 

280.05 

-212.41 

99.40 

5194.75 

-- 

16 

J 

5572.08 

-544.36 

280.05 

-212.41 

99.40 

5194.75 

-- 

16 

6 

5213.68 

-386.10 

154.57 

-96.98 

33.49 

4918.66 

— 

16 

5 

5213.68 

-386.10 

154.57 

-96.98 

33.49 

4918.66 

-- 

16 

4 

4888.72 

-272.77 

78.09 

-57.13 

6.23 

4663.14 

-- 

16 

3 

4888.72 

-272.77 

78.09 

-37.13 

6.23 

4663.14 

-- 

16 

2 

4598.22 

-196.35 

34.98 

-8.93 

-2.76 

4425.15 

— 

16 

1 

4598.09 

-196.31 

34.96 

-8.92 

-2.77 

4425.05 

-- 

16 

0 

4344.76 

-150.33 

13.88 

2.25 

-4.18 

4206.39 

-- 

16 

-1 

4342.94 

-149.43 

13.53 

2.42 

-4.26 

4205.20 

-- 

16 

-2 

4137.73 

-132.16 

7.65 

4.73 

-3.02 

4014.93 

-- 

16 

-3 

4122.07 

-123.83 

4.86 

5.67 

-3.46 

4005.30 

-- 

16 

-4 

3992.13 

-134.13 

9.35 

3.82 

-1.72 

3869.44 

-- 

16 

-5 

3924.18 

-107.31 

1.52 

5.79 

-2.51 

3821.67 

-- 

16 

-6 

3874.30 

-124.67 

7.53 

3.78 

-1.52 

3759.41 

-- 

i 


Tabl 


J 

T 

er 

*4 

*6 

16 

-7 

3725.71 

-88.54 

-0.92 

16 

-8 

3716.63 

-95.31 

1.15 

16 

-9 

3505.84 

-66.01 

-2.62 

16 

-10 

3505.01 

-6/. 00 

-2.36 

10  (Continued) 


*8 

*10 

et 

Observed 

4.92 

-1.83 

3639.34 

3657.1 

4.31 

-1.63 

3625.16 

3640.3 

3.46 

-1.20 

3439.47 

3439.7 

3.39 

-1.18 

3437.86 

3437.7 

From  the  work  in  chapter  3,  expectations  are  that  due  to  the 


strong  coupling  between  the  vibrational  and  rotational  motions,  the 
various  coefficients  in  equation  (5.3)  will  exhibit  a strong  dependence 
on  whichever  set  of  states  is  used  for  fitting  purposes.  In  these  cal- 
culations, the  fit  is  to  energy  levels  instead  of  spectral  lines.  Al- 
though there  are  far  more  spectral  lines  than  levels,  both  sets  of  data 
contain  exactly  the  same  information  on  the  physical  structure  of  the 
system.  Owing  to  the  large  number  of  lines,  the  standard  deviations 
obtained  in  a least  squares  analysis  will  be  much  smaller  than  that 
obtained  from  fits  to  energy  levels.  However,  since  the  additional 
equations  can  all  be  obtained  from  those  equations  used  to  fit  to 
energies,  it  is  unlikely  that  the  value  of  the  rotation  constants,  etc., 
would  change  significantly.  Furthermore,  the  coefficients  derived 
from  spectral  lines  will  be  much  more  strongly  correlated  than  those 
derived  from  an  energy  level  analysis  and,  consequently,  less  reliable. 

In  table  11,  the  coefficients  that  arise  fTom  a least  squares 
analysis  of  high  t states  are  compared  to  those  derived  from  low  t 
states.  More  precisely,  the  following  68  levels  were  used  for  the  high 
t fit:  5 1 — 5 5 j 61-65;  7 1 - 77 ; 81-83;  9 1 — 9 9 ; 10i-10iq  and  llj-lln.  The 
low  t fit  used  l-i»lo»li»  2_ 2 , 2_ 1 ; 3. 3 - 3_ i ; 4_4-4_i;  5_ 5- 5_ 1 ; 6-6-6-1; 

7_ 7- 7_  1 ; 8_ q- 8_  1 ; 9_g-9-i;  10_io-10-i  and  ll_n-ll_i.  In  both  cases, 
these  coefficients  yielded  energy  eigenvalues  that  deviated  by  only  a 
few  hundredths  of  a wave  number  from  the  measured  energy  levels  used 
for  fitting. 

An  examination  of  table  11  reveals  that  most  of  the  fourth- 
order  and  virtually  all  of  the  higher  order  distortion  coefficients  of 
the  low  and  high  t fits  deviate  from  one  another  by  more  than  a standard 
deviation.  For  example,  at  the  center  of  their  distributions,  the  high 
and  low  x differ  by  a factor  of  2.  For  &j,  the  difference  is  an 
order  of  magnitude  and  for  5^,  two  orders  of  magnitude.  Also,  there  is 
a sign  difference.  These  deviations  are  serious  because  the  fourth- 
order  energy  is  generally  large  and  implies  serious  errors  would  be 
present  for  any  wave  function  generated  from  these  constants.  In  general, 
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the  higher  order  coefficients  deviate  from  one  to  four  orders  of  magnitude 
with  sign  changes.  The  value  of  these  coefficients  is  then  completely 
determined  by  the  immediate  levels  used  for  fitting  and  do  not  reflect 
the  physical  structure  of  the  molecule. 

The  cause  of  these  difficulties  lies  in  the  fact  that  the  rotation- 
vibration  energy  of  the  high  r states  is  much  greater  than  that  of  the 
low  t levels.  Since  no  effort  has  been  made  to  remove  this  energy  from 
the  distortion  coefficients,  they  strongly  reflect  its  presence.  A 
second  difficulty  lies  in  the  fact  that  the  separation  between  the  ro- 
tational levels  (i.e.,  between  different  t's  of  a fixed  J and  a different 
vibrational  band)  is  now  comparable  to  the  vibrational  energy  splitting. 

For  example,  the  vibrational  energy  difference  between  the  ground  and 
first  excited  bending  mode  is  about  1600  cm-1 . For  the  rotational 
splitting 

EV  = l,  J = 9,  t = 9 " EV  = 0,  J = 9,  t = 5 

- 1600  cm  1 “ 800  cm"1 , 

which  is  half  the  vibrational  energy  difference.  These  coefficients 
depend  strongly  on  the  energy  difference  between  the  rotation-vibration 
levels  that  are  coupled  together  by  the  nonrigid  terms  of  the  Hamiltonian. 
Consequently,  they  will  depend  strongly  on  the  structure  of  the  rotational 
state  when  the  rotational  energy  splitting  becomes  comparable  to  the 
vibrational  spacings. 

In  table  12,  a comparison  is  made  between  the  high  and  low  J 
levels.  In  particular,  a fit  was  made  to  70  levels  in  the  J = 1 to  8 
states  and  to  all  of  the  J = 10,  11,  12  states.  Deviations  again  arise 
in  the  fourth  order,  as  in  the  case  of  high  and  low  t.  However,  owing  to 
the  large  standard  deviations,  less  detail  can  be  obtained  from  this  table. 

PREDICTABILITY 

The  ability  of  equation  (5.3)  to  predict  the  positions  of  new 
energy  levels  has  been  examined.  Specifically  a comparison  was  made  of  the 
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predicted  rotational  levels  that  arise  from  the  Watson  Hamiltonian  [i.e., 
equation  (5.3)]  for  the  following  sets  of  rotation-distortion  coefficients: 
(1)  low  x;  (2)  high  t;  (3)  low  J ; (4)  high  J,  (5)  reference  12. 

Our  results  show  that  serious  discrepancies  appear  in  all  cases, 
although  reference  12  is  far  superior  to  the  others. 

Low  r 

First,  the  results  of  fitting  to  the  low  x states  will  be  discussed. 
In  table  13  are  listed  the  predicted  values  versus  the  observed  values 
for  the  different  rotational  states  through  J = 16.  Note  that 

those  levels  that  were  chosen  for  fitting  purposes  are  reproduced  to 
better  than  0.05  cm'1.  First,  consider  the  high  x states  for  J <.  11. 

(Recall  that  the  low  x states  for  J < 11  were  used  for  fitting  pur- 

poses.) An  examination  of  table  13  reveals  that  deviations  on  the  order 
of  a few  tenths  of  a wave  number  have  already  appeared  by  the  J = 3 
and  4 states.  For  the  52  and  53  states,  deviations  of  2 cm'1  occur;  and 
for  the  54  and  5s  states,  there  is  a discrepancy  between  theory  and 
experiment  of  about  10  cm  1 . Similarly,  for  the  63,  64  states,  deviations 
of  about  0.7  cm  1 appear;  and  for  the  65  and  6g  rotational  levels,  errors 
on  the  order  of  40  cm  1 have  appeared.  This  trend  continues,  i.e., 
greater  and  greater  discrepancies  between  theory  and  experiment  arise 
in  the  J = 8 and  J = 9 rotational  multiplets.  Finally,  by  J = 10, 

there  is  no  resemblance  between  the  predicted  energies  and  the  observed 
ones.  For  example,  the  10j  and  IO2  states  deviate  by  9 cm'1;  for  IO7 
and  IO3,  the  predicted  energies  are  negative,  a situation  that  arose  for 
the  case  of  a rotating  Morse  oscillator. 

Next,  consider  the  low  x levels  for  J > 12.  An  examination  of 

table  13  reveals  that  12_i2  - 12.2  rotational  energies  are  in  reasonable 
agreement  with  experiment.  However,  beyond  this,  the  situation  rapidly 
deteriorates.  Similarly,  the  states  13_i 3 - 13. 3 are. fairly  well  pre- 
dicted, although  the  13.2  level  deviates  by  15  cm'1  from  the  experiment. 

By  J = 14,  serious  descrepancies  arise  by  x = -6,  for  J = 15  by 
x = -7,  and  for  J = 16  by  x = -14. 
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TEST  FOR  UNIQUENESS  OF  LEAST  SQUARES  FITTING 
FOR  HIGH  AND  LOW  t 


Low  t High  r 


A 

2.7861234 

X 

10 1 

± 

2.1908565 

X 

10'2 

2.790118 

X 

IO1 

t 

1.1406407 

X 

10-2 

B 

1.45217S1 

X 

10* 1 

♦ 

3.8477284 

X 

10*  3 

1.4445913 

X 

IO1 

+ 

8.1378196 

X 

10-2 

C 

9.2816871 

± 

3. 1S8S214 

X 

IO*3 

9.3003991 

1.0102989 

X 

IO-1 

V 

1 . 27S1778 

X 

10*  3 

3.9871204 

X 

10*  5 

-2.1136767 

X 

IO*3 

2.7358925 

X 

IO-3 

V* 

-4.2498711 

X 

10"  3 

1.2848932 

X 

10* 3 

-4.6654608 

X 

10*  4 

8.7506809 

X 

io- 3 

1.63807S0 

X 

io-2 

± 

6.9675567 

X 

io*3 

3.4500498 

X 

10'2 

5.6726837 

X 

io- 3 

SK 

4.9866468 

X 

IO'4 

± 

2.3341136 

X 

10*5 

4.9825552 

X 

10- 3 

± 

2.7474778 

X 

10*  3 

SK 

1.S4228S3 

X 

io*4 

1.3868588 

X 

io*3 

-5.4694406 

X 

IO"2 

+ 

3.5702885 

X 

io*2 

BJ 

6.6118469 

X 

10*7 

2.2555794 

X 

10*7 

1.5868118 

X 

IO-4 

4; 

7.8205148 

X 

IO'5 

H _ 

2.4312591 

X 

10*5 

9.3611777 

X 

10*6 

-2.1034987 

X 

IO*3 

1.0772050 

X 

io*3 

JK 

SKJ 

2.1972752 

X 

10*4 

4- 

1.5221769 

X 

io*4 

4.9642308 

X 

10*  3 

■¥ 

2.5955161 

X 

10*  3 

-2.6745956 

X 

io*3 

+ 

1.0754196 

X 

10*  3 

-2.9120505 

X 

10* 3 

■¥ 

1.5904144 

X 

io- 3 

h. 

1.7076044 

X 

10- 7 

•f 

1.2714262 

X 

10*7 

-7.6465990 

X 

IO*5 

± 

3.7979067 

X 

io*5 

0 

h 

7. S151086 

X 

IO*5 

± 

3.7290733 

X 

10*6 

1.5419958 

X 

10* 3 

* 

7.7448027 

X 

10*4 

JK 

hK 

-2.3733418 

X 

10*  4 

2.5585506 

X 

io*4 

-S. 4848512 

X 

10* 3 

± 

2.9108518 

X 

io*3 

ljk 

-2.1781492 

X 

io*6 

+ 

7.5353344 

X 

io*7 

4.2607054 

X 

10*s 

2.5026737 

X 

10*6 

L„  „ 

1.4896433 

X 

10*  5 

± 

9.1730935 

X 

io*6 

-1.1822368 

X 

IO*5 

± 

6.9779067 

X 

io*5 

KJ 

lk 

7.1554451 

X 

10*5 

± 

5.6518559 

X 

10*  5 

7.2549200 

X 

IO*6 

i 

4.4860100 

X 

10*6 

1.7274188 

X 

10*5 

± 

1.3524660 

X 

10*  5 

8.8590460 

X 

10*s 

5.4370631 

X 

io*6 

-1.4597506 

X 

10*  5 

7.5929448 

X 

IO*7 

5.1953952 

X 

10- : 10 

± 

1.7756653 

X 

io*10 

-2.3539468 

X 

IO*7 

± 

2.2110777 

X 

IO*7 

4.5974880 

X 

IO*9 

+ 

2.8869232 

X 

10*9 

► 
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Table  12 

TEST  FOR  UNIQUENESS  OF  LEAST  SQUARES  FITTING 
FOR  J = 1-8,  10,  11,  12 


J * 10. 

11.  12 

« J - 

1- 

8 

A 

2.8021076 

X 

10 1 

± 

1.2584273 

X 

10" 1 

2.7877399 

X 

IO1 

4- 

5.7776983 

X 

10- 3 

B 

1 . 447S6S4 

X 

101 

4 

3.8022468 

X 

IO'2 

1.4521195 

X 

IO1 

± 

3.5283291 

X 

10- 3 

C 

9.2749972 

♦ 

1.9043611 

X 

IO'2 

9.2772397 

4 

3.1810502 

X 

10- 3 

V 

9 . 40496S3 

X 

10-“ 

± 

2.4951000 

X 

10-u 

1.2206764 

X 

10" 3 

♦ 

8.5179333 

X 

10- 5 

V* 

-4.4687892 

X 

10*3 

4 

2.0059901 

X 

io- 3 

-6.0803482 

X 

io- 3 

± 

5.4899506 

X 

io*4 

3.  S623288 

X 

10‘2 

± 

2.1873472 

X 

io'3 

3.2391491 

X 

10'2 

± 

9.0682603 

X 

io-4 

3.5035443 

X 

IO*4 

4 

1.8738260 

X 

IO"4 

5.1759149 

X 

10‘4 

4 

7.4186159 

X 

io- 5 

5k 

1.6797911 

X 

io- 3 

4 

3.3470801 

X 

io*4 

1.2575220 

X 

io- 3 

4 

1.2739147 

X 

10' 3 

HJ 

-3.8357118 

X 

io- 7 

± 

9.2518603 

X 

io-7 

5.7229858 

X 

IO"7 

4 

1.2563280 

X 

IO"6 

bjk 

-2.5619922 

X 

10*6 

± 

7.4562113 

X 

10"6 

-1.4591014 

X 

10*  5 

t 

1.6255732 

X 

10"  5 

BKJ 

1.2740388 

X 

10*5 

4 

2.0188034 

X 

10' 5 

-2.7046471 

X 

10"  5 

4 

4.0917037 

X 

io-5 

h 

1.3581069 

X 

IO"4 

4 

2.4762050 

X 

10' 5 

1.4235874 

X 

10"4 

- 

5.3898739 

X 

io-5 

hJ 

-2.7959194 

X 

10"7 

6.9202722 

X 

10'7 

3.0319408 

X 

10'7 

4 

6.7186041 

X 

io- 7 

hJK 

1.8733722 

X 

io-7 

± 

3.5403482 

X 

io-6 

1.2472925 

X 

io- 5 

4 

1.9690785 

X 

io"5 

hK 

4.3887342 

X 

io-5 

± 

2.2086231 

X 

io-5 

-6.8265611 

X 

10' 5 

4 

2.8141018 

X 

10‘4 

ljk 

-3.5071943 

X 

10*  7 

2.5076621 

X 

io*7 

3.1812650 

X 

io- 6 

4 

4.4365160 

X 

10- 6 

9.2501065 

X 

io*7 

± 

7.5724814 

X 

io*7 

-8.4370406 

X 

io-° 

4 

1.2577229 

X 

10- 5 

KJ 

lk 

-1.1867769 

X 

io- 6 

± 

S. 0924892 

X 

io*7 

4.9919688 

X 

io*6 

4 

8.2803728 

X 

IO*6 

lK 

-1.0525886 

X 

io-6 

+ 

7.1830927 

X 

10"  7 

4.2204533 

X 

10'6 

4 

1.2441923 

X 

io-s 

PK 

1.1888523 

X 

io-9 

t 

5.1584160 

X 

io-10 

-3.0137686 

X 

10* 10 

t 

4.8525514 

X 

IO*9 

Px 

6.4596572 

X 

io-10 

± 

9.0045005 

X 

io- 10 

3.7693584 

X 

10"8 

4.1545549 

X 

10- 3 
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Table  13 

PREDICTABILITY  OF  VARIOUS  FITS 


J_ 

T 

Obs. 

Low  t 

High  t 

Low  J 

High  J 

1 

1 

23.791 

23.80 

23.76 

23.80 

23.79 

1 

0 

37.137 

37.13 

37.08 

37.13 

37.02 

1 

+1 

42.362 

42.37 

42.39 

42.37 

42.26 

2 

-2 

70.088 

70.10 

70.03 

70.09 

70.09 

2 

-1 

79.496 

79.49 

79.41 

79.49 

79.40 

2 

0 

95.174 

95.17 

95.16 

95.18 

95.08 

2 

1 

134.902 

134.92 

134.86 

134.89 

134.48 

2 

2 

136.1611 

136.18 

136.19 

136.16 

135.75 

3 

-3 

136.765 

136.78 

137.03 

136.76 

136.74 

3 

-2 

142.274 

142.28 

142.62 

142.28 

142.20 

3 

-1 

173.370 

173.37 

173.21 

173.38 

173.29 

3 

0 

206.303 

206.33 

206.15 

206.30 

205.96 

3 

1 

212.162 

212.18 

212.18 

212.16 

211.85 

3 

2 

285.217 

284.90 

285.21 

285.21 

284.37 

3 

3 

285.421 

285.10 

285.42 

285.41 

284.67 

4 

-4 

222.050 

222.08 

224.06 

222.06 

222.01 

4 

-3 

224.844 

224.84 

227.05 

224.84 

224.77 

4 

-2 

275.494 

275.50 

275.09 

275.51 

275.42 

4 

-1 

300.367 

300.39 

300.19 

300.37 

300.11 

4 

0 

315.777 

315.81 

315.68 

315.79 

315.57 

4 

1 

382.521 

382.36 

382.43 

382.52 

381.89 

4 

2 

383.837 

383.69 

383.84 

383.84 

383.23 

4 

3 

488.110 

485.55 

488.12 

488.12 

486.84 

4 

4 

488.126 

485.58 

488. 15 

488.14 

486.87 

5 

-5 

325.399 

325.38 

332.70 

325.36 

325.29 

5 

-4 

326.620 

326.62 

334.25 

326.64 

326.56 

5 

-3 

399.50 

399.47 

399.01 

399.47 

399 . 37 

5 

-2 

416.202 

426.21 

416.49 

416.23 

416.04 

5 

-1 

446.56 

446.53 

445.93 

446.53 

446.38 

5 

0 

593.961 

593.92 

503.78 

503.99 

503.56 

5 

1 

508.86 

508.78 

508.83 

508.82 

508.44 

5 

0 

U 

610. 100 

608.51 

610.09 

610.12 

609.24 

5 

3 

610.38 

608.75 

610.32 

610.35 

609.47 

S 

4 

742.10 

731.51 

742.10 

742.10 

740.54 

5 

5 

742.11 

731.52 

742.10 

742.11 

740.55 

r 


Table  13 

(Continued) 

£ 

T 

Obs. 

Low  t 

High  t 

Low  J 

High  J 

6 

-6 

446.704 

446.72 

466.89 

446.73 

446.65 

6 

-5 

447.30 

447.25 

467.71 

447.29 

447.20 

6 

-4 

542.890 

432.94 

543.75 

542.91 

542.82 

6 

-3 

552.96 

552.91 

555.06 

552.93 

552.81 

6 

-2 

602.762 

602.79 

601.22 

602.79 

602.68 

6 

-1 

649.02 

648.99 

648.67 

649.01 

648.77 

6 

0 

661 . 534 

661.57 

661.33 

661.56 

661.37 

6 

1 

756.76 

755.91 

756.68 

756.75 

756.24 

6 

2 

757.78 

756.99 

757.78 

757.80 

757.30 

6 

3 

888.65 

881.75 

888.61 

888.61 

887.62 

6 

4 

888.59 

881.79 

888 . 6 j 

888.64 

887.66 

6 

5 

104S.11 

1004.05 

1045.08 

1045.07 

1043.59 

6 

6 

1045.03 

1004.05 

1045.09 

1045.07 

1043.59 

7 

n 
— / 

586.28 

586.26 

633.36 

586.29 

586.21 

7 

-6 

586.43 

586.48 

633.81 

586.53 

586.44 

7 

-5 

704.20 

704.27 

716.11 

704.21 

704.17 

7 

-4 

709.50 

709.62 

717.18 

709.60 

709.58 

7 

-3 

782.41 

782.42 

779.53 

782.44 

782.33 

7 

-2 

816.65 

816.72 

816.40 

816.73 

816.65 

7 

-1 

842.36 

842.39 

841.10 

842.34 

842.29 

7 

0 

927.76 

927.49 

927.59 

927.78 

927.59 

7 

1 

931.23 

931.03 

931.27 

931.26 

931.09 

7 

2 

1059.68 

1055.97 

1059.67 

1059.77 

1059.19 

7 

3 

1059.89 

1056.18 

1059.83 

1059.96 

1059.38 

7 

4 

1216.39 

1186.53 

1216.24 

1216.20 

1215.33 

7 

5 

1216.39 

1186.54 

1216.24 

1216.21 

1215.34 

7 

6 

1394.85 

1220.1-6 

1394.85 

1394.80 

1393.83 

/ 

7 

1394.85 

1220.16 

1394.85 

1394.80 

1393.83 

8 

-8 

744.09 

744.08 

842.60 

744.11 

744.05 

8 

-7 

744.14 

744.17 

842.85 

744.20 

744.14 

8 

-6 

882.97 

882.98 

901.78 

882.85 

882.92 

8 

'5 

885.64 

885.64 

906.19 

885.54 

885.64 

8 

-4 

983.09 

982.96 

979.27 

983.00 

982.89 

8 

-3 

1006.14 

1006.16 

1006.68 

1006.14 

1006.18 

8 

-2 

1050.20 

1050.18 

1046.58 

1050.14 

1050.14 

8 

-1 

1122.78 

1122.76 

1122.14 

1122.72 

1122.80 

8 

0 

1131.88 

1131.89 

1131.50 

1131.77 

1131.84 

8 

1 

1255.19 

1253.87 

1255.21 

1255.49 

1255.14 

8 

2 

1255.98 

1254.65 

1255.95 

1256.25 

1255.89 

8 

3,4 

1411.59 

1391.69 

1411.66 

1412.28 

1411.30 

8 

5 , 6 

1590.70 

1445.84 

1590.81 

1590.78 

1590.19 

k 
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Table  13  (Continued) 


J_ 

T 

Obs . 

Low  T 

High  t 

Low  J 

High  J 

8 

7,8 

1789.09 

1068.80 

1789.11 

1789.76 

1788.90 

9 

-9 

920.20 

920.19 

1110.24 

920.16 

920.18 

9 

-8 

920.22 

920.23 

1110.38 

920.20 

920.22 

9 

-7 

1079.16 

1079.20 

1125.76 

1078.97 

1079.19 

9 

-6 

1080.51 

1080.49 

1128.51 

1080.23 

1080.49 

9 

-5 

1202.04 

1202.04 

1200.51 

1202.15 

1202.00 

9 

-4 

1216.37 

1216.31 

1220.33 

1216.27 

1216.37 

9 

-3 

1283.02 

1282.96 

1275.59 

1282.99 

1282.93 

9 

-2 

1341.03 

1341.06 

1339.36 

1340.89 

1341.12 

9 

-1 

1360.56 

1360.45 

1358.43 

1360.19 

1360.37 

9 

0 

1475.14 

1475.11 

1474.83 

1475.54 

1475.25 

9 

1 

1477.46 

1477.51 

1477.38 

1477.91 

1477.55 

9 

2 

1631.44 

1618.68 

1631.33 

1633.23 

1631.36 

9 

3 

1631.58 

1618.85 

1631.14 

1633.39 

1631.50 

9 

4,5 

1810.76 

1691.12 

1810.70 

1813.20 

1810.49 

9 

6,7 

2009.99 

1080.49 

2010.05 

2010.57 

2009.79 

9 

8,9 

2225.56 

-402.17 

2225.6 

2231.26 

2226.16 

10 

-10 

1114.59 

1114.59 

1458.71 

1114.40 

1114.59 

10 

-9 

1114.59 

1114.60 

1458.79 

1114.42 

1114.60 

10 

-8 

1293.22 

1293.17 

1393.15 

1292.81 

1293. 18 

10 

-7 

1293.80 

1293.79 

1394.85 

1293.39 

1293.78 

10 

-6 

1438.19 

1438.15 

1446.20 

1438.49 

1438.17 

10 

-5 

1446.23 

1446.28 

1460.28 

1446.34 

1446.29 

10 

-4 

1538.31 

1538.28 

1526.10 

1538.54 

1538.26 

10 

-3 

1581.53 

1581.54 

1578.48 

1581.35 

1581.60 

10 

-2 

1616.49 

1616.6 

1610.79 

1616.51 

1616.57 

10 

-1 

1719.36 

1719.35 

1717.66 

1719.57 

1719.11 

10 

0 

1724.80 

1725.42 

1724.5 

1725.67 

1725.05 

10 

1 

1875.24 

1866.80 

1875.06 

1878.93 

1875.34 

10 

2 

1875.72 

1867.38 

1875.53 

1879.49 

1875.83 

10 

3,4 

2054.55 

1953.92 

2054.44 

2062.19 

2054.55 

10 

5,6 

2254.55 

1663.37 

2254.54 

2262.43 

2254.39 

10 

7,8 

2471.59 

-4.64 

2471.67 

2475.64 

2471.59 

10 

9,10 

2702.09 

-5678.57 

2702.21 

2728.54 

2702.83 

11 

-10 

1327.25 

1327.23 

1919.04 

1326.76 

1327.24 

11 

-9 

1525.02 

1525.02 

1720.42 

1524.58 

1525.02 

11 

-8 

1525.31 

1525.33 

1721.46 

1524.83 

1525.29 

11 

-7 

1690.85 

1690.84 

1722.93 

1691.81 

1690.92 

11 

-6 

1695.24 

1695.27 

1732.74 

1695.85 

1695.21 

11 

-5 

1813.47 

1813.51 

1797.41 

1814.47 

1813.54 

Table  13  (Continued) 


£ 

X 

Obs . 

Low  x 

High  x 

Low  J 

High  J 

11 

-4 

1843.32 

1843.26 

1839.63 

1843.52 

1843.24 

11 

-3 

1899.21 

1899.17 

1886.34 

1899.61 

1899.19 

11 

-2 

1968.08 

1968.11 

1982.30 

1987.22 

1986.14 

11 

-1 

1999.34 

1999.31 

1997.22 

2000.58 

1999.26 

11 

0 

2143.01 

2135.5 

2142.30 

2149.28 

2143.04 

11 

1 

2144.46 

2137.18 

2144.11 

2150.90 

2144.48 

11 

2 

2322.20 

2232.40 

2321.9 

2337.89 

2322.18 

11 

3 

2322.25 

2232.59 

2322.93 

2338.80 

2322.27 

11 

4,5 

2522.46 

1975.18 

2522.45 

2536.28 

2522.50 

11 

6,7 

2740.73 

391.39 

2740.90 

2762.64 

2740.64 

11 

8,9 

2973.07 

-5100.52 

2973.43 

2991.54 

2973.13 

11 

10,11 

3216.6 

-20510.2 

3216.74 

3309.34 

3216.59 

12 

-11 

1558.07 

1558.10 

2632.95 

1557.10 

1558.07 

12 

-10 

1774.75 

1774.75 

2131.17 

1774.45 

1774.75 

12 

-9 

1774.88 

1774.91 

2131.77 

1774.56 

1774.86 

12 

-8 

1960.38 

1960.25 

2041.81 

1962.71 

1960.40 

12 

-7 

1962.60 

1962.66 

2048.45 

1964.66 

1962.55 

12 

-6 

2106.7 

2106.30 

2091.60 

2109.04 

2166.43 

12 

-5 

2124.84 

2125.26 

2125.03 

2126.98 

2125.05 

12 

-4 

2205.95 

2205.80 

2182.65 

2207.84 

2206.12 

12 

-3 

2275.65 

2274.84 

2267.23 

2278.29 

2275.60 

12 

-2 

2300.94 

2299.92 

2294.21 

2303.62 

2300.89 

12 

-1 

2434.14 

2424.51 

2431.84 

2444.38 

2434.17 

12 

0 

2437.84 

2428.58 

2437.32 

2448.32 

2437.84 

12 

1 

2613.26 

2525.07 

2612.72 

2640.55 

2613.22 

12 

2 

2613.49 

2525.69 

2613.04 

2640.98 

2613.52 

12 

3,4 

2813.94 

2290.85 

2813.58 

2863.07 

2813.94 

12 

5,6 

3033.17 

776.45 

3033.08 

3095.21 

3033.15 

12 

7,8 

3267.2 

-4552.03 

3267.49 

3322.66 

3267.23 

12 

9,10 

3512.8 

-19625.6 

3513.24 

3577.09 

3512.78 

12 

11,12 

3767.1 

-57005.3 

3767.35 

4039.92 

3767.11 

13 

-12 

1806.94 

1807.12 

3355.04 

1805.35 

1807.07 

13 

-11 

2042.5 

2042.29 

2657.93 

2042.62 

2042.32 

13 

-10 

2042.5 

2042.38 

2658.27 

2042.67 

2042.37 

13 

-9 

2247.0 

2246.61 

2419.79 

2252.17 

2246.86 

13 

-8 

2248.24 

2247.97 

2424.12 

2253.06 

2247.89 

13 

-7 

2415.95 

2415.15 

2515.60 

2421.83 

2415.30 

13 

-6 

2426.0 

2426.62 

2440.29 

2431.94 

2426.14 

13 

-5 

2534.14 

2533.83 

2498.00 

2539.84 

2534.87 

13 

-4 

2586.5 

2584.84 

2571.60 

2593.10 

2586.63 

13 

-3 

2629.54 

2626.94 

2613.01 

2635.38 

2629.86 

45 


Table  13  (Continued) 


J_ 

T 

Obs . 

Low  T 

High  t 

Low  J 

High  J 

13 

_2 

2748.4 

2733.45 

2741.74 

2764,80 

2748.86 

13 

-1 

2756.61 

2742.40 

2755.00 

2772.98 

2756.64 

13 

0 

2927.38 

2830.87 

2926.10 

2970.87 

2927.54 

13 

1 

2928.45 

2832.61 

2927.60 

2971.94 

2928.39 

13 

2,3 

3128.25 

2606.58 

3127.55 

3214.14 

3128.63 

13 

4,5 

3348.2 

1141.36 

3347.83 

3476.44 

3348.99 

13 

6,7 

3584.0 

-4051.36 

3583.93 

3729.26 

3585.09 

14 

-13 

2073.81 

2074.26 

4455.38 

2071.36 

2074.18 

14 

-12,-11 

2328.2 

2327.55 

3344.54 

2329.42 

2327.71 

14 

-10 

2551.0 

2550.00 

2881.32 

2561.44 

2550.48 

14 

-9 

2551.5 

2550.81 

2884.04 

2561.83 

2550.95 

14 

-3 

2740.5 

2739.75 

2777.00 

2754.14 

2739.99 

14 

-7 

2745.5 

2746.62 

2795.93 

2759.37 

2745.77 

14 

-6 

2883.5 

2880.40 

2833.20 

2895.31 

2882.81 

14 

-5 

2919.5 

2915.40 

2896.48 

2932.91 

2918.34 

14 

-4 

2983.6 

2978.52 

2950.65 

2996.26 

2984.97 

14 

-3 

3085.0 

3062.08 

2069.50 

3112.04 

2085.14 

14 

-2 

2101.65 

3079.68 

3095.90 

3126.69 

2101.92 

14 

-1 

3264.2 

3149.09 

3260.45 

3330. 3S 

3265.01 

14 

0 

3266.36 

3153.50 

3265.55 

3332.42 

3267.16 

14 

1 

3465.18 

2919.25 

3463.73 

3601.09 

3466.33 

14 

2 

3465.4 

2920.39 

3464.02 

3601.39 

2466.51 

15 

-14 

2358. S8 

2359.48 

5922.20 

2355.02 

2359.39 

15 

-13,-12 

2631.6 

2630.24 

4247. SS 

2635.20 

2630.78 

15 

-11 

2872.56 

2870.39 

3460.38 

2892.11 

2871.31 

15 

-10 

2872.9 

2870.89 

3462.01 

2892.28 

2871.52 

15 

-9 

3081.2 

3080.57 

3195.62 

3108.79 

3080.44 

15 

-8 

3084.2 

3084.83 

3208.91 

3111.28 

3083.37 

15 

-7 

3252.0 

3243.24 

3193.34 

3275.73 

3247.56 

15 

-6 

3277.0 

3265.96 

3246.43 

3300.36 

3269.89 

15 

-5 

3365.0 

3351.38 

3304.39 

3387.88 

3364.21 

15 

-4 

3446.0 

3410.08 

3412.86 

3489.23 

5444.01 

15 

-3 

3473.0 

3440.70 

3456.91 

3512.13 

3474.67 

15 

-2 

3624.0 

3479.27 

3613 . 14 

3722.01 

3625.56 

15 

-1 

3628.7 

3489.49 

3626.89 

3724.61 

3630.45 

15 

0 

3824.8 

3226.47 

3820.93 

4027.40 

3827.28 

15 

1 

3862.1 

3229.82 

3822.42 

4027.86 

3827.75 

15 

2,3 

4045.8 

1778.31 

4042.73 

4400.08 

4050.58 

16 

-16,-15 

2661.2 

2662.73 

7864.79 

2656.20 

2662.68 

16 

-14,-13 

2953.0 

2950.21 

5440.64 

2960.64 

2951.52 

46 
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Table  13  (Continued) 


J_ 

T 

Obs . 

Low  T 

High  t 

Low  J 

High  J 

16 

-12,-11 

3211.5 

3207.66 

4202.60 

3246.26 

2662.68 

16 

-10 

3437.7 

3438.43 

3695.47 

3489.79 

3437.12 

16 

-9 

3439.7 

3441.21 

3704.34 

3490.89 

3438.54 

16 

-8 

3640.3 

3621.81 

3588.72 

3685.03 

3627.63 

16 

-7 

3657.1 

3636.46 

3631.32 

3699.97 

3640.54 

y 


High  t 

Neat,  examine  the  high  t fit.  An  examination  of  table  13 
shows  that  deviations  of  a few  wave  numbers  arise  for  low  J levels, 
such  as  4_4,  4_3,  5.5,  5.4.  For  6-g,  6.5,  it  is  ^ 20  cm-1,  and  this 
trend  continues  (getting  worse)  for  all  the  low  t states,  e.g.,  for 
10_io, it  is  over  300  wave  numbers.  These  difficulties  arise  from  and 
reflect  the  strong  rotation-vibration  correlations  that  exist  for  high 
x levels,  which  in  turn  manifest  themselves  in  the  various  coefficients. 

Low  J 

The  Hamiltonian  was  chosen  so  that  the  J =1  through  8 states 
were  fit  to  'v  0.05  cm-1.  An  examination  of  table  13  reveals  that  dif- 
ficulties manifest  themselves  almost  immediately.  For  example,  in  the 
92  state,  there  is  a deviation  of  2 cm-1;  and  by  9g>  it  is  6 cm-*.  The 
situation  continues  to  worsen,  e.g.,  10}  is  off  by  2 cm’1,  lOjg  by  26  cm-1. 
For  J = 1,  difficulties  arise  by  11-2  (1  cm-1);  and  by  lln,  it  is 

100  cm'1.  Thus,  the  ability  of  this  fit  to  predict  high  states  is 
virtually  nonexistent. 

High  J 

The  J - 10,  11,  and  12  levels  were  fit.  With  few  ex- 

ceptions, all  of  the  states  were  fit  to  better  than  0.05  cm-1.  The 
pattern  is  very  similar  to  the  situation  that  arose  for  the  rotating 
Morse  oscillator  when  the  levels  J = 10  to  20  were  fit;  namely,  the 

low-lying  levels  are  reasonably  well  predicted.  However,  there  was 
immediate  difficulty  for  the  higher  states,  i.e.,  136  deviated  by  1 cm-1, 
143  by  3 cm'1,  15s  by  3 cm"1,  153  by  5 cm'1,  16. 8 by  13  cm'1  and  I6-7  by 
17  cm-1 . 

Microwave  spectral  Data 

This  work  (ref.  1)  is  a spectroscopic  fit  to  the  spectral  data  that 
arise  from  microwave  transitions  between  the  low  energy  rotational  states 
through  J - 10.  The  results  of  predictions  for  this  case  are  quali- 

tatively similar,  although  quantitatively  superior  to  the  predictions 
of  the  other  fits.  An  examination  of  tables  9 and  10  reveals  the  fol- 
lowing features.  All  of  the  rotational  levels  through  89  are  predicted 
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to  within  a few  hundredths  of  a wave  number.  Beyond  this,  more 
serious  difficulties  arise.  For  example,  for  the  92  through  9g  levels, 
deviations  on  the  order  of  a few  tenths  of  a wave  number  arise  between 
the  predicted  levels  and  the  measured  ones.  Furthermore,  some  of  the 
high  t states  of  the  J =10  multiplet  are  off  by  almost  half  of  a 
wave  number.  For  the  remaining  levels,  a pattern  that  is  similar  to 
that  for  high  and  low  t,  as  well  as  for  high  and  low  J appears,  al- 
though this  case  is  quantitatively  more  accurate.  Within  the  J = 11 

series,  deviations  on  the  order  of  several  wave  numbers  arise;  and  for 
the  12n  and  12x2  states  (where  S’io  > 100  cm-1),  discrepancies  of  over 
30  cm-1  appear.  At  the  moment,  enough  data  are  not  available  beyond 
the  J - 12  states  to  thoroughly  examine  further  predictions  of  the 

microwave  fit.  However,  it  is  undoubtedly  safe  to  assume  that  de- 
viations continue  to  increase  and,  in  addition,  that  they  appear  for 
lower  t numbers  as  J increases. 


p 


■■ 


SECTION  VI 

DISCUSSION  AND  CONCLUSIONS 

In  this  paper  we  have  examined  the  usefulness  and  reliability  of 
power  series  expansion  techniques  (as  embodied  in  Watson's  rotational 
Hamiltonian)  in  elucidating  the  rotational  structure  of  light  asym- 
metric highly  anharmonic  molecules,  particularly  water.  This  question 
has  been  examined  within  the  framework  of  both  a model  system  (i.e.,  the 
rotating  Morse  oscillator)  as  well  as  in  the  water  molecule  itself.  In 
particular,  for  water,  the  Hamiltonian  has  terms  of  order  P10  and  involves 
18  perturbation  coefficients  in  addition  to  the  three  rotational  constants. 
In  both  cases,  the  expansion  techniques  proved  useful  only  for  the  low- 
lying  rotational  levels.  For  such  states,  the  rotation-vibration  inter- 
action is  relatively  weak,  and  a power  series  expansion  based  on  the  pic- 
ture of  a slightly  deformed  rotor  is  valid.  This  is  borne  out  by  the  rapid 
convergence  of  the  various  perturbation  terms  displayed  in  table  9.  Similar 
behavior  was  exhibited  by  the  rotating  Morse  oscillator.  On  the  other 
hand,  in  attempts  to  apply  the  same  procedure  to  the  high- lying  rotational 
states,  serious  difficulties  immediately  appeared.  For  those  states,  the 
molecular  rotation  velocity  is  large,  and  both  the  centripetal  and 
coriolis  (for  polyatomic  molecules)  forces  are  large  and,  in  turn,  deform 
the  molecule  in  a significant  manner.  Since  these  deformations  represent 
huge  perturbations  away  from  the  rigid  rotor  structure,  many  rotation- 
distortion  terms  are  required  to  correctly  describe  the  physical  state 
of  the  molecule.  To  date,  calculations  are  limited  to  tenth-order  per- 
turbations (i.e.,  P10  terms),  and  the  convergence  results  exhibited  in 
tables  9 and  10  show  that  this  is  not  enough,  i.e.,  P12,  P14,  etc.,  type 
terms  are  required.  The  type  of  difficulties  that  arise  depends  on  the 
data  used  to  determine  the  various  rotation  and  rotation-distortion  co- 
efficients. Thus,  data  that  arise  from  the  low-lying  levels  cannot  be  used 
to  accurately  predict  the  level  structure  for  the  higher  states.  In  addi- 
tion, although  the  first  few  rotation-distortion  coefficients  are  accurately 
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determined,  the  higher-order  ones  may  not  be.  This  is  due  to  the  fact 
that  the  low-lying  levels  are  only  weakly  deformed  from  the  rigid  rotor 
structure  and,  consequently,  are  not  a sensitive  measure  of  P8  and  P10 
type  terms.  Evidence  of  this  behavior  appears  in  tables  10  and  13. 
Attempts  to  circumvent  these  problems  by  fitting  to  higher  rotational 
levels  results  in  a new  host  of  difficulties.  In  particular,  the 
rotational  constants  deviate  considerably  from  their  expected  values, 
and  this  gives  rise  to  an  incorrect  picture  of  the  geometric  structure 
of  the  molecules.  In  addition,  it  also  appears  that  one  cannot  predict 
the  rotation  energies  of  the  low-lying  levels  (see  table  13)  for  HjO, 
which  implies  that  all  of  the  lower-order  rotation  distortion  coef- 
ficients are  in  serious  error.  This  is  clearly  brought  out  in  tables 
11  and  12.  In  the  case  of  the  rotating  Morse  oscillator,  the  low-lying 
level  structures  were  predicted  reasonably  accurately,  although  the 
various  rotation  and  rotation-distortion  coefficients  were  in  strong 
disagreement  with  the  low  energy  fit.  The  physical  basis  for  these 
difficulties  lies  in  the  fact  that  in  attempts  to  fit  the  distortion 

. 

constants  to  the  high-lying  levels,  higher-order  rotation-vibration 
correlations  (i.e.,  12th  etc.)  are  folded  into  the  lower  order  coef- 
ficients. As  a consequence,  the  rotational  wave  functions  will  bear 

it 

little  resemblance  to  reality. 

Another  difficulty  that  arises  is  due  to  the  large  number  of 
perturbation  coefficients  and  the  limited  amount  of  convergent  levels 
to  fit  these  coefficients.  This  causes  many  of  the  coefficients  that 
are  determined  in  a least  squares  manner  to  be  highly  correlated.  In 
turn,  serious  doubts  are  cast  as  to  whether  or  not  these  numbers 
really  represent  the  physical  quantities  that  they  should.  A related 
problem  arises  in  trying  to  include  higher  order  distortion  effects, 
i.e.,  P12,  P14,  ...  terms  in  the  Watson  Hamiltonian.  In  particular, 
there  are  many  more  unknown  coefficients  to  determine  and  only  a limited 
number  of  energy  levels.  Although  there  are  far  more  spectral  lines 
than  energy  levels,  these  lines  contain  no  more  information  on  the 
structure  of  the  molecule  than  do  the  rotational  energies.  Hence, 


although  there  are  more  data  points  to  fit,  they  are  highly  correlated; 
this,  in  turn,  gives  rise  to  strong  correlations  among  the  higher 
order  coefficients  thus  reducing  their  reliability. 

Another  feature  that  should  be  pointed  out  is  that  the  rotation 
constants  depend  on  the  energy  denominators  (ST.T  - T ),  which  can 

VtJ  j V e/T  i 

be  replaced  by  ( Eyj  - Ey,j  ) only  for  the  low-lying  states.  For  higher 
energy  levels,  the  rotational  splitting  can  be  on  the  order  of  the  vibra- 
tional spacing.  The  rotation  and  rotation-distortion  constants  then 
depend  sensitively  on  the  particular  rotational  states  used  for  fitting 
purposes,  and  these  will  not  be  unique,  i.e.,  they  will  vary  with  the 
set  of  levels  used  as  input  data.  This  is  clearly  displayed  in  tables 
7,  11,  and  12. 

Therefore,  although  this  approach  will  work  for  low  energy 
rotational  levels,  it  does  not  work  and,  on  the  theoretical  grounds  just 
mentioned,  cannot  work  for  higher  J or  high  r states.  If  reliable  wave 
functions  need  to  be  obtained  for  the  higher  rotational  states  (J  z 10) 
in  addition  to  the  lower  ones,  another  approach  would  be  required. 
Finally,  in  order  to  obtain  the  rotational  level  structure  in  a different 
vibrational  band,  the  entire  process  would  have  to  be  repeated  since 
the  various  rotation  and  rotation-distortion  coefficients  will  vary 
significantly  with  the  vibrational  state  of  the  molecule. 


APPENDIX 

THE  EFFECTS  OF  MATRIX  DIAGONALIZATION  ON  THE  DISTORTION  ENERGIES 


In  this  appendix  we  examine  the  effects  of  matrix  diagonalization 
on  the  various  rotational  distortion  energies,  i.e.,  E 4,  Eg>  etc. 
convenience,  decompose  the  rotational  Hamiltonian  into  a rigid  rotor 
piece  JK'o  plus  a distortion  term  .V  where 

+ <*8  + JK'lO  (Al) 

In  equation  (Al) , JtC’g,  ...  refer  to  the  fourth,  sixth,  ...  order 

distortion  terms. 

Note  that  the  results  of  table  9 show  that  effects  of  matrix 
diagonalization  on  the  total  energy  are  small,  which  justifies  using 
perturbation  theory  in  the  analysis  given  below.  The  leading  correction 
to  the  rotational  eigenvector  due  to  3?'  is 
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In  equation  (A2) , |^T>^  D refers  to  the  matrix  diagonalized  eigenvector 
and  |cl  > to  the  rigid-rotor  eigenvector.  The  effects  of  matrix  diagonal- 
ization, for  any  distortion  term,  say  the  fourth, 
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To  estimate  the  size  of  these  terms,  note  that  the  perturbing  Hamiltonian 
has  all  the  symmetries  consequently,  it  will  connect  the  state  J t 


only  to  those  states  J , such  that 
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t'  = t ± 4nj  n = 1,2,.,.,  (A4) 

Now,  due  to  the  large  spacings  between  those  different  states  that 
are  connected  by  1 it  follows  that  the  major  contributions  (A2)  and 
(A3)  arise  from  the  n = 1 term.  As  a specific  example,  we  consider 
the  1010  state.  Then,  the  dominant  contribution  arises  from  the  106 
state  and  since  from  the  definition  of 


= y &E 
L n 

n = 4 


so  that  for  1 Oj q , 
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from  which  it  follows  that 

| < 1010  | JIT*  | 106  > | * 5 cm’1 

To  evaluate  the  leading  contributions  to  6£’(4),...  in  general, 
expect 

(A7) 
(AS) 

and  the  various  corrections  are 
I | < 2.76  cm-1 
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